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Preface 



Of the five editions of this text, this is the third edition that I have prepared. In 
doing so, I have attempted to adhere to the late Jerry Marion’s original purpose of 
producing a modem and reasonably complete account of the classical mechanics 
of particles, systems of particles, and rigid bodies for physics students at the ad- 
vanced undergraduate level. The purpose of the book continues to be threefold: 

1. To present a modem treatment of classical mechanical systems in such a way 
that the transition to the quantum theory of physics can be made with the 
least possible difficulty. 

2. To acquaint the student with new mathematical techniques wherever possi- 
ble, and to give him/her sufficient practice in solving problems so that the 
student may become reasonably proficient in their use. 

3. To impart to the student, at the crucial period in the student’s career be- 
tween “introductory” and “advanced” physics, some degree of sophistication 
in handling both the formalism of the theory and the operational technique 
of problem solving. 

After a firm foundation in vector methods is presented in Chapter 1, further 
mathematical methods are developed in the textbook as the occasion demands. 
It is advisable for students to continue studying advanced mathematics in sepa- 
rate courses. Mathematical rigor must be learned and appreciated by students of 
physics, but where the continuity of the physics might be disturbed by insisting 
on complete generality and mathematical rigor, the physics has been given 
precedence. 



Changes for the Fifth Edition 

The comments and suggestions of many users of Classical Dynamics have been in- 
corporated into this fifth edition. Without the feedback of the many instructors 
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Special Feature 

The author has kept one popular feature of Jerry Marion’s original book: the ad- 
dition of historical footnotes spread throughout. Several users have indicated 
how valuable these historical comments have been. The history of physics has 
been almost eliminated from present-day curricula, and as a result, the student is 
frequently unaware of the background of a particular topic. These footnotes are 
intended to whet the appetite and to encourage the student to inquire into the 
history of his field. 



Teaching Aids 

Teaching aids to accompany the textbook are available online at 
http://info.brookscole.com/thornton . The Instructor’s Manual (ISBN 0-534- 
40898-2) contains solutions to all the end-of-chapter problems in addition to 
Transparency Masters of selected key figures from the text. This password- 
protected resource is easily printable in .pdf format. To receive your password, 
just go to the above website and register; a username and password will be sent to 
you once the information you have provided is verified. The verification 
procedure ensures that you are an instructor teaching this course. If you are not 
able to download the Instructor’s Manual files and would like a printed copy sent 
to you, please contact your local sales representative. If you do not know who 
your sales representative is, please visit www.brookscole.com, and click on the 
Find your Rep tab, which is located at the top of the web page. Please do not dis- 
tribute the Instructor’s Manual to students, or post the solutions on the Internet. 
Students are not permitted to access the Instructor’s Manual. 



Student Solutions Manual 

A Student Solutions Manual by Stephen T. Thornton, which contains solutions to 
25% of the problems, is available for sale to the students. Instructors are encour- 
aged to order the Student Solutions Manual for their students to purchase at the 
school bookstore. To package the Student Solutions Manual with the text, use 
ISBN 0-534-08378-1, or to order the Student Solutions Manual separately use ISBN 
0-534-40897-4. Students can also purchase the manual online at the publisher’s 
website www.brookscole.com/physics. 
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PREFACE 



who have used this text, it would not be possible to produce a textbook of signif- 
icant value to the physics community. After the extensive revision for the fourth 
edition, the changes in this edition have been relatively minor. Only a few re- 
arrangements of material have been made. But several examples, especially nu- 
merical ones, and many end-of-chapter problems have been added. Users have 
not wanted extensive changes in the topics covered, but more examples for stu- 
dents and a wider range of problems are always requested. 

A strong effort continues to be made to correct the problem solutions avail- 
able in the Instructor and Student Solutions Manuals. I thank the many users 
who sent comments concerning various problem solutions, and many of their 
names are listed below. Answers to even-numbered problems have again been in- 
cluded at the end of the book, and the selected references and general biblio- 
graphy have been updated. 



Course Suitability 

The book is suitable for either a one-semester or two-semester upper level (jun- 
ior or senior) undergraduate course in classical mechanics taken after an intro- 
ductory calculus-based physics course. At the University of Virginia we teach a 
one-semester course based mostly on the first 12 chapters with several omissions 
of certain sections according to the Instructor’s wishes. Sections that can be 
omitted without losing continuity are denoted as optional, but the instructor can 
also choose to skip other sections (or entire chapters) as desired. For example, 
Chapter 4 (Nonlinear Oscillations and Chaos) might be skipped in its entirety 
for a one-semester course. Some instructors choose not to cover the calculus of 
variations material in Chapter 6. Other instructors may want to begin with 
Chapter 2, skip the mathematical introduction of Chapter 1, and introduce the 
mathematics as needed. This technique of dealing with the mathematics intro- 
duction is perfectly acceptable, and the community is divided on this issue with a 
slight preference for the method used here. The textbook is also suitable for a 
full academic year course with an emphasis on mathematical and numerical 
methods as desired by the instructor. 

The textbook is appropriate for those who choose to teach in the traditional 
manner without computer calculations. However, more and more instructors 
and students are both familiar and adept with numerical calculations, and much 
can be learned by doing calculations where parameters can be varied and real- 
world conditions like friction and air resistance can be included. I decided be- 
fore the 4th edition to leave the choice of method to the instructor and/or stu- 
dent to choose the computer techniques to be used. That decision has been 
confirmed, because there are many excellent software programs (including 
Mathematica, Maple, and Mathcad to mention three) available to use. In addi- 
tion, some Instructors have students write computer programs, which is an im- 
portant skill to obtain. 
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CHAPTER 




Matrices, Vectors, 
and Vector Calculus 



1.1 Introduction 

Physical phenomena can be discussed concisely and elegantly through the use of 
vector methods.* In applying physical “laws” to particular situations, the results 
must be independent of whether we choose a rectangular or bipolar cylindrical 
coordinate system. The results must also be independent of the exact choice of 
origin for the coordinates. The use of vectors gives us this independence. A 
given physical law will still be correctly represented no matter which coordinate 
system we decide is most convenient to describe a particular problem. Also, the 
use of vector notation provides an extremely compact method of expressing 
even the most complicated results. 

In elementary treatments of vectors, the discussion may start with the state- 
ment that “a vector is a quantity that can be represented as a directed line seg- 
ment.” To be sure, this type of development will yield correct results, and it is 
even beneficial to impart a certain feeling for the physical nature of a vector We 
assume that the reader is familiar with this type of development, but we forego 
the approach here because we wish to emphasize the relationship that a vector 
bears to a coordinate transformation. Therefore, we introduce matrices and ma- 
trix notation to describe not only the transformation but the vector as well. We 
also introduce a type of notation that is readily adapted to the use of tensors, al- 
though we do not encounter these objects until the normal course of events re- 
quires their use (see Chapter 11). 



*Josiah Willard Gibbs (1839-1903) deserves much of the credit for developing vector analysis 
around 1880-1882- Much of the present-day vector notation was originated by Oliver Heaviside 
(1850-1925), an English electrical engineer, and dates from about 1893. 
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1 / MATRICES, VECTORS, AND VECTOR CALCULUS 



We do not attempt a complete exposition of vector methods; instead, we 
consider only those topics necessary for a study of mechanical systems. Thus in 
this chapter, we treat the fundamentals of matrix and vector algebra and vector 
calculus. 



1.2 Concept of a Scalar 

Consider the array of particles shown in Figure 1-la. Each particle of the array is 
labeled according to its mass, say, in grams. The coordinate axes are shown so 
that we can specify a particular particle by a pair of numbers (.*, y). The mass M 
of the particle at (x, y) can be expressed as M(x, y ); thus the mass of the particle 
at x = 2, y = 3 can be written as M (x = 2, y = 3) — 4. Now consider the axes ro- 
tated and displaced in the manner shown in Figure 1-lb. The 4 g mass is now lo- 
cated at x' = 4, y' = 3.5; that is, the mass is specified by M (x r = 4, y’ = 3.5) = 4. 
And, in general, 

M(x, y) = M{x',y') (1.1) 

because the mass of any particle is not affected by a change in the coordinate 
axes. Quantities that are invariant under coordinate transformation — those that obey 
an equation of this type — are termed scalars. 

Although we can describe the mass of a particle (or the temperature, or the 
speed, etc.) relative to any coordinate system by the same number, some physical 
properties associated with the particle (such as the direction of motion of the 
particle or the direction of a force that may act on the particle) cannot be speci- 
fied in such a simple manner. The description of these more complicated quan- 
tities requires the use of vectors. Just as a scalar is defined as a quantity that re- 
mains invariant under a coordinate transformation, a vector may also be defined 
in terms of transformation properties. We begin by considering how the coordi- 
nates of a point change when the coordinate system rotates around its origin. 




FIGURE 1-1 An array of particles in two different coordinate systems. 














1 .3 COORDINATE TRANSFORMATIONS 
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1.3 Coordinate Transformations 

Consider a point Pwith coordinates (x 1; x 2 , x s ) with respect to a certain coordi- 
nate system.* Next consider a different coordinate system, one that can be gen- 
erated from the original system by a simple rotadon; let the coordinates of the 
point Pwith respect to the new coordinate system be (x(, x 2 , X 3 ). The situation is 
illustrated for a two-dimensional case in Figure 1-2. 

The new coordinate x[ is the sum of the projection of Xj onto the x[-axis 
(the line Oa) plus the projection of x 2 onto the x(-axis (the line ab + be); that is, 

(1.2a) 

Od — de, but the 

x 2 = — x t sin 6 + x 2 cos 6 
( it \ 

= Xjcosl — + 6 I + x 2 cos 6 (1.2b) 

Let us introduce the following notation: we write the angle between the 
xj-axis and the X]-axis as (x{, xj), and in general, the angle between the x'-axis 
and the x 7 -axis is denoted by (x-, xj). Furthermore, we define a set of numbers 
A„ by 

A,y = cos(x', Xj) (1.3) 



x{ = Xj cos 6 + x 2 sin 6 



rr 



= Xj cos 6 + x 2 cos^— — dj 

The coordinate x 2 is the sum of similar projections: x 2 
line de is also equal to the line Of. Therefore 



* 2 -axis 




FIGURE 1-2 The position of a point Fean be represented in two coordinate systems, 
one rotated from the other. 



*We label axes as atj, x. 2 , instead of x, y, z to simplify the notation when summations are performed. 
For the moment, the discussion is limited to Cartesian (or rectangular) coordinate systems. 
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1 / MATRICES, VECTORS, AND VECTOR CALCULUS 



Therefore, for Figure 1-2, we have 

An = cos(x(, *i) = cos 9 



A 12 = cos(xJ, x 2 ) = cosl — — 9 J = sin 9 



A 2 i = cos(x 2 , »i) = cos^— + 6 J = -sin 9 

A 22 = cos(x 2 , x 2 ) = cos 9 / 

The equations of transformation (Equation 1.2) now become 

x( = Xi cos(x{, x^ + x 2 cos(x{, x 2 ) 
AnXi f A 12 x 2 

x 2 = XiCOS(x 2 , Xj) + x 2 cos(x 2 , x 2 ) 

~ A 2J Xi T A 22 x 2 

Thus, in general, for three dimensions we have 

x{ = AnXi + Ai 2 x 2 + A 13 x 3 ^ 
x 2 = A 2 iXi "H A 22 x 2 "I - A 23 x 3 / 

X3 = A 3 iXi + A S2 x 2 + A33X3 J 
or, in summation notation. 



3 

x' = 2 a ijXj, i- 1, 2, 3 

j— i 1 J 



The inverse transformation is 



(1.4) 

(1.5a) 

(1.5b) 

( 1 . 6 ) 

(1.7) 



Xi = x{ COS(x(, Xi) + X 2 COS(x 2 , X 2 ) + X3 COS (x 3 , Xi) 

AjjXj f A 2 iX 2 -f A31X3 

or, in general, 

( 1 . 8 ) 

The quantity A,y is called the direction cosine of the x'-axis relative to the 
x ; -axis. It is convenient to arrange the A y into a square array called a matrix. The 
boldface symbol A denotes the totality of the individual elements \ tJ when 
arranged as follows: 





Once we find the direction cosines relating the two sets of coordinate axes, 
Equations 1.7 and 1.8 give the general rules for specifying the coordinates of a 
point in either system. 

When A is defined this way and when it specifies the transformation proper- 
ties of the coordinates of a point, it is called a transformation matrix or a rota- 
tion matrix. 
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EXAMPLE l.l 



A point Pis represented in the (xj, x 2 , x 3 ) system by P( 2, 1, 3). In another coor- 
dinate system, the same point is represented as P(x{, x 2 , x 3 ) where x 2 has been 
rotated toward x 3 around the Xj-axis by an angle of 30° (Figure 1-3). Find the 
rotation matrix and determine P(x{, x 2 , x 3 ). 




FIGURE 1-3 Example 1.1. A point Pis represented in two coordinate -systems, one 
rotated from the other by 30°. 



Solution. The direction cosines A t y can be determined from Figure 1-3 using 
the definition of Equation 1 .3. 



A 12 

^13 

A 2 i 

a 22 

^23 

^31 

^32 

^33 



COS(x(, Xi) 
cos(x(, x 2 ) 
COS (Xj, x 3 ) 
COS(x 2> Xj) 
cos(x 2 , x 2 ) 
cos ( x 2 , x 3 ) 
cos ( x 3 , Xj) 
cos(x 3 , x 2 ) 
cos(x 3 , x 3 ) 



A = 



cos(0°) - 1 
cos (90°) = 0 
cos (90°) = 0 
cos (90°) = 0 
cos (30°) = 0.866 
cos(90° - 30°) = cos(60°) = 0.5 
cos(90°) = 0 
cos (90° + 30°) = -0.5 
cos (30°) = 0.866 
/l 0 0 

0 0.866 0.5 

\0 -0.5 0.866, 



and using Equation 1.7, P(x(, x 2 , x 3 ) is 



x{ = A n x! + A 12 x 2 4- A 13 x 3 — Xi = 2 

x 2 = A 2 jXj + A 22 x 2 + A 23 x 3 = 0.866x 2 + 0.5x 3 = 2.37 

x 3 = A 3 }X} + A 32 x 2 + A 33 x 3 = — 0.5x 2 4- 0.866x s = 2.10 
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Notice that the rotation operator preserves the length of the position vector, 
r = Vx? + xf + x 3 = A/xJ 2 + x 2 2 + x 3 2 = 3.74 



1.4 Properties of Rotation Matrices* 

To begin the discussion of rotation matrices, we must recall two trigonometric 
results. Consider, as in Figure l-4a, a line segment extending in a certain direc- 
tion in space. We choose an origin for our coordinate system that lies at some 
point on the line. The line then makes certain definite angles with each of the 
coordinate axes; we let the angles made with the x r , x 2 -, x 3 -axes be a, /3, y. The 
quantities of interest are the cosines of these angles; cos a, cos /3, cos y. These 
quantities are called the direction cosines of the line. The first result we need is 
the identity (see Problem 1-2) 

cos 2 a + cos 2 /3 + cos 2 y = 1 (1*10) 

Second, if we have two lines with direction cosines cos a, cos (3, cos y and cos a', 
cos /S', cos y ', then the cosine of the angle 6 between these lines (see Figure l-4b) 
is given (see Problem 1-2) by 

cos 6 = cos a cos a' + cos /8 cos /S' + cos y cos y' (1*11) 

With a set of axes x lt x 2 , x 3 , let us now perform an arbitrary rotation about 
some axis through the origin. In the new position, we label the axes x{, x 2 , x 3 . 




(a) (b) 



FIGURE 1^4 (a) A line segment is defined by angles (a, (3, y) from the coordinate axes. 

(b) Another line segment is added that is defined by angles ( a (3', y'). 



*Much of Sections 1.4-1 .13 deals with matrix methods and transformation properties and will not be 
needed by the reader until Chapter 1 1 . Hence the reader may skip these sections until then if de- 
sired. Those relations absolutely needed — scalar and vector products, for example — should already 
be familiar from introductory courses. 
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The coordinate rotation may be specified by giving the cosines of all the angles 
between the various axes, in other words, by the Ay. 

Not all of the nine quantities Ay are independent; in fact, six relations exist 
among the Ay, so only three are independent. We find these six relations by 
using the trigonometric results stated in Equations 1.10 and 1.11. 

First, the x{-axis may be considered alone to be a line in the (x 1? x 2 , x 3 ) coor- 
dinate system; the direction cosines of this line are (A u , A 12 , A 13 ). Similarly, the di- 
rection cosines of the x 2 -axis in the (aq, x 2 , x 3 ) system are given by (A 21 , A 22 , A 23 ). 
Because the angle between the xj-axis and the x 2 -axis is tt/2, we have, from 
Equation 1.11, 

A h A 2] T Aj 2 A 22 T A 13 A 23 — cos 0 = cos (tt/2) = 0 



or* 

A [ j A 2 j = 0 

And, in general, 

Xa*A*,-=0, i±k (1.12a) 

i 

Equation 1.12a gives three (one for each value of i or ft) of the six relations 
among the Ay. 

Because the sum of the squares of the direction cosines of a line equals unity 
(Equation 1.10), we have for the x(-axis in the (xj, x 2 , x 3 ) system, * 



^11 + ^12 + ^13 — 1 



or 



and, in general, 



/ ) 



^)j ^-kj 1) t A 



(1.12b) 



which are the remaining three relations among the A y . 

We may combine the results given by Equations 1.12a and 1.12b as 




(M3) 



where 8 tk is the Kronecker delta symboT 



. _ r°, 

** u. 



if i ^ ft 
if i ~ ft 



(1-14) 



The validity of Equation 1.13 depends on the coordinate axes in each of the 
systems being mutually perpendicular. Such systems are said to be orthogonal, 



*A11 summations here are understood to run from 1 to 3. 
^Introduced by Leopold Kronecker (1823—1891). 
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and Equation 1.13 is the orthogonality condition. The transformation matrix A 
specifying the rotation of any orthogonal coordinate system must then obey 
Equation 1.13. 

If we were to consider the x r axes as lines in the x\ coordinate system and 
perform a calculation analogous to our preceding calculations, we would find 
the relation 



2a,vA,^ Sj k 



(1.15) 



The two orthogonality relations we have derived (Equations 1.13 and 1.15) 
appear to be different. (Note: In Equation 1.13 the summation is over the second 
indices of the A,-,, whereas in Equation 1.15 the summation is over the first in- 
dices.) Thus, it seems that we have an overdetermined system: twelve equations 
in nine unknowns.* Such is not the case, however, because Equations 1.13 and 
1.15 are not actually different. In fact, the validity of either of these equations 
implies the validity of the other. This is clear on physical grounds (because the 
transformations between the two coordinate systems in either direcdon are 
equivalent), and we omit a formal proof. We regard either Equation 1.13 or 1.15 
as providing the orthogonality relations for our systems of coordinates. 

In the preceding discussion regarding the transformation of coordinates 
and the properties of rotation matrices, we considered the point P to be fixed 
and allowed the coordinate axes to be rotated. This interpretation is not unique; 
we could equally well have maintained the axes fixed and allowed the point to 
rotate (always keeping constant the distance to the origin) . In either event, the 
transformation matrix is the same. For example, consider the two cases illustrated 
in Figures l-5a and b. In Figure l-5a, the axes X] and x 2 are reference axes, and 
the x(- and x 2 -axes have been obtained by a rotation through an angle 6. 



*2 





*1 

(b) 



FIGURE 1-5 (a) The coordinate axes X],x 2 are rotated by angle 6, but the point P 

remains fixed, (b) In this case, the coordinates of point Pare rotated 
to a new point P', but not the coordinate system. 



♦Recall that each of the orthogonality relations represents six equations. 
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Therefore, the coordinates of the point P with respect to the rotated axes may 
be found (see Equations 1.2a and 1.2b) from 

x[ = Xj cos 9 + x 2 sin 9 

*2 = ~~ sin 9 + x 2 cos 0 

However, if the axes are fixed and the point Pis allowed to rotate (as in Figure 
l-5b) through an angle 9 about the origin (but in the opposite sense from that 
of the rotated axes), then the coordinates of P' are exactly those given by 
Equation 1.16. Therefore, we may elect to say either that the transformation acts 
on the point giving a new state of the point expressed with respect to a fixed co- 
ordinate system (Figure l-5b) or that the transformation acts on the frame of ref- 
erence (the coordinate system), as in Figure l-5a. Mathematically, the interpreta- 
tions are entirely equivalent. 



(1.16) 



1.5 Matrix Operations* 



The matrix A given in Equation 1.9 has equal numbers of rows and columns and 
is therefore called a square matrix. A matrix need not be square. In fact, the co- 
ordinates of a point may be written as a column matrix 



x = 



( Xj\ 

x 2 

\*s/ 



(1.17a) 



or as a row matrix 



X = (*! *2 *s) 



(1.17b) 



We must now establish rules to multiply two matrices. These rules must be 
consistent with Equations 1.7 and 1.8 when we choose to express the x t and the 
x[ in matrix form. Let us take a column matrix for the coordinates; then we have 
the following equivalent expressions: 




x' 


= 2. 
j 


KjXj 




(1.18a) 


x' 


= Ax 






(1.18b) 


/A n 


a 12 


Al 3 ^ 


/ 




( A 2 i 


a 22 


A23 


x 2 


(1.18c) 


\A 31 


A32 


A33/ 


\* 3 / 





X] — AjjXj T Aj 2 ^2 A" A13X3 I 

X 2 — A 21 X l "P A 2 2#2 ”P A23X3 / 

x 3 = AjiXj ^■ 32 x 2 A33X3 / 



(1.18d) 



*The theory of matrices was first extensively developed by A. Cayley in 1855, but many of these ideas 
were the work of Sir William Rowan Hamilton (1805-1865), who had discussed “linear vector opera- 
tors” in 1852. The term matrix was first used byj. J. Sylvester in 1850. 
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Equations 1.18a-d completely specify the operation of matrix multiplication 
for a matrix of three rows and three columns operating on a matrix of three 
rows and one column. (To be consistent with standard matrix convention we 
choose X and x' to be column matrices; multiplication of the type shown in 
Equation 1.18c is not defined if x and x' are row matrices.)* We must now ex- 
tend our definition of multiplication to include matrices with arbitrary numbers 
of rows and columns. 

The multiplication of a matrix A and a matrix B is defined only if the num- 
ber of columns of A is equal to the number of rows of B. (The number of rows of 
A and the number of columns of B are each arbitrary.) Therefore, in analogy 
with Equation 1.18a, the product AB is given by 



C = AB 

0/ — t AB] y = 2 A ik B k j 



As an example, let the two matrices A and B be 



( 1 . 19 ) 




-2 

-3 

b 




* f 



h j 



We multiply the two matrices by 



AB 




-2 

-3 




b c 

e f 



h j 



(1.20) 



The product of the two matrices, C, is 

C = AB = h~ 2d+2g M- 2e+2h 3c-2/+2A a21) 

\4a — 3d + 5g 4b — 3e + 5 h 4c — 3 /+ 5jJ 

To obtain the C t] element in the rth row andjth column, we first set the two 
matrices adjacent as we did in Equation 1.20 in the order A and then B. We then 
multiply the individual elements in the «th row of A, one by one from left to 
right, times the corresponding elements in the ^th column of B, one by one 
from top to bottom. We add all these products, and the sum is the Qj element. 
Now it is easier to see why a matrix A with m rows and n columns must be multi- 
plied times another matrix B with n rows and any number of columns, say p. The 
result is a matrix C of m rows and p columns. 



“Although whenever we operate on X with the A matrix the coordinate matrix X must be expressed 
as a column matrix, we may also write x as a row matrix (x,, x 2 , x 3 ), for other applications. 
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Solution. We follow the example of Equations 1.20 and 1.21 to multiply the two 
matrices together. 




/-2 + 1+9 -4 + 2 + 12\ / 8 10\ 

AB = 2 + 2 + 12 4 + 4 + 16 = 16 24 

\ 1-3- 12 2-6-16/ V — 14 -20/ 



The result of multiplying a 3 X 3 matrix times a 3 X 2 matrix is a 3 X 2 matrix. 



It should be evident from Equation 1.19 that matrix multiplication is not 
commutative. Thus, if A and B are both square matrices, then the sums 

B kj and ^LB lk A k , 

are both defined, but, in general, they will not be equal. 



EXAM PI. Iv 1.3 



Show that the multiplication of the matrices A and B in this example is non- 
commutative. 

Solution. If A and B are the matrices 






then 
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but 



thus 




AB A BA 



1 . 6 Further Definitions 



A transposed matrix is a matrix derived from an original matrix by interchange 
of rows and columns. We denote the transpose of a matrix A by A*. According to 
the definition, we have 



Evidendy, 



4 = 






(A*)* = A 



( 1 . 22 ) 

(1.23) 



Equation 1.8 may therefore be written as any of the following equivalent expres- 
sions: 



X.' 

j 

X; ij Xj 

} 

x = A'x' 



*1^1 




fhl 


A 2 i 


Asi^ 


*2 


= 


Aig 


a 22 


A32 


xj 




\^1S 


^23 


^33^ 



(x[\ 

x 2 

W 



(1.24a) 

(1.24b) 

(1.24c) 

(1.24d) 



The identity matrix is that matrix which, when multiplied by another matrix, 
leaves the latter unaffected. Thus 



that is, 



1A = A, B1 = B 






(1.25) 



Let us consider the orthogonal rotauon matrix A for the case of two dimensions: 
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A n 

^12 



A21 j 

a 22 / 

AnA 2 ] + A 12 A 22 

A|] + A99 



Then 

AA' = f^* A ’ 2 ' 

\A 21 a 22 , 

_( Afi 4- Af 2 

\A 2 iA n + A 22 A 12 a 2 ] t a 22 

Using the orthogonality relation (Equation 1.13), we find 

Afi Af 2 — A 21 + A| 2 = 1 

A 21 A 11 + A 22 A 12 = A]]A2i “J" A]2A22 = 0 
so that for the special case of the orthogonal rotation matrix A we have* 

'1 0 



AA 



i — 



.0 1. 



= 1 



(1.26) 



The inverse of a matrix is defined as that matrix which, when multiplied by 
the original matrix, produces the identity matrix. The inverse of the matrix A is 
denoted by A -1 : 



A A 1 = 1 



(1.27) 



By comparing Equations 1.26 and 1.27, we find 




for orthogonal matrices 



(1.28) 



Therefore, the transpose and the inverse of the rotation matrix A are identical. 
In fact, the transpose of any orthogonal matrix is equal to its inverse. 

To summarize some of the rules of matrix algebra: 



1 . Matrix multiplication is not commutative in general: 

AB / BA (1.29a) 

The special case of the multiplication of a matrix and its inverse is commu- 
tative: 



AA 1 = A 1 A = 1 (1.29b) 

The identity matrix always commutes: 

1A = A1 = A (1.29c) 

2. Matrix multiplication is associative: 

[AB]C = A[BC] (1.30) 

3. Matrix addition is performed by adding corresponding elements of the two 
matrices. The components of C from the addition C = A + B are 

q, = Aj + q (i.3i) 

Addition is defined only if A and B have the same dimensions. 



*This result is not valid for matrices in general. It is true only for orthogonal matrices. 
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1.7 Geometrical Significance 
of Transformation Matrices 

Consider coordinate axes rotated counterclockwise* through an angle of 90° 
about the x 3 -axis, as in Figure 1-6. In such a rotation, xj = x 2 , x 2 = — x lt x 3 = x 3 . 
The only nonvanishing cosines are 

cos(x{, x 2 ) = 1 = A 12 

COs(x 2 , Xj) = — 1 = A 2 i 
cos (x 3 , x 3 ) = 1 = A 33 

so the A matrix for this case is 

( 0 1 °\ 

A, = I — 1 0 0 

V 0 0 1/ 

Next consider the counterclockwise rotation through 90° about the x r axis, 
as in Figure 1-7. We have x\ — x l5 x 2 — x 3 , x 3 — — x 2 , and the transformation ma- 
trix is 

(l 0 0\ 

A 2 = 0 0 1 

\0 -1 0/ 



To find the transformation matrix for the combined transformation for rota- 
tion about the x 3 -axis, followed by rotation about the new xj-axis (see Figure 
1-8), we have 



and 



or 



/ 



*1 




x' = AjX 

x" = A 2 x' 



X = A 2 AjX 



/I 


0 


°) 


/ 0 


1 


0\ 


1 x x ) 


0 


0 


1 


-1 


0 


0 


x 2 


\0 


-1 


0/ 


\ 0 


0 


l) 


\xj 



(° 


1 


0\ 


^x/ 


° 


0 


1 


x 2 


\l 


0 


0/ 


\x 3 / 



(1.32a) 

(1.32b) 

(1.33a) 




(1.33b) 



*We determine the sense of the rotation by looking along the positive portion of the axis of rotation 
at the plane being rotated. This definition is then consistent with the “right-hand rule,” in which the 
positive direction of advance of a right-hand screw when turned in the same sense. 




1 .7 GEOMETRICAL SIGNIFICANCE OF TRANSFORMATION MATRICES 



15 




FIGURE 1-6 Coordinate system x lt x%, x s is rotated 90° counter-clockwise (ccw) 
about the x 3 -axis. This is consistent with the right-hand rule of 
rotation. 




x, x{ 

FIGURE 1-7 Coordinate system x 1; x 2 , *3 is rotated 90° ccw about the x r axis. 



*3 




FIGURE 1-8 Coordinate system x lt x 2 , x 3 is rotated 90° ccw about the x 3 -axis followed 
by a 90° rotation about the intermediate x,'-axis. 
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Therefore, the two rotations already described may be represented by a single 
transformation matrix: 





(0 


l 


0\ 


A s — AgA] — 


0 


0 


1 




U 


0 


0/ 



(1.34) 



and the final orientation is specified by x" — x 2 , x 2 = x 5 , x'i, = x x . Note that the 
order in which the transformation matrices operate on X is important because 
the multiplication is not commutative. In the other order, 





0 


1 


°\ 


f l 


0 


°\ 




-1 


0 


0 


0 


0 


1 




0 


0 


1/ 


v0 - 


1 


0/ 


( 


0 


0 


1 N 


V 








-1 


0 


0 


1 ^ A3 






\ 


0 


-1 


0, 


/ 







(1.35) 



and an entirely different orientation results. Figure 1-9 illustrates the different 
final orientations of a parallelepiped that undergoes rotations corresponding to 
two rotation matrices A a , X b when successive rotations are made in different 
order. The upper portion of the figure represents the matrix product X B X A , and 
the lower portion represents the product X A X B . 




FIGURE 1-9 A parallelepiped undergoes two successive rotations in different order. 
The results are different. 
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Next, consider the coordinate rotation pictured in Figure 1-10 (which is the 
same as that in Figure 1-2). The elements of the transformation matrix in two 
dimensions are given by the following cosines: 



COs(xJ, Xi) 
cos(x^, x 2 ) 

COS (x 2 , Xi) 
cos(x 2 , x 2 ) 

Therefore, the matrix is 



= cos 6 = A n 




= cos 6 = A 2 2 



( cos 6 sin 6 

— sin 6 cos 0 



(1.36a) 



If this rotation were a three-dimensional rotation with x 3 = x 3 , we would 
have the following additional cosines: 



cos(x(, x 3 ) = 0 = A 13 
cos(x 2 , x 3 ) 0 A 23 

cos(x 3 , x 3 ) = 1 = A 33 
cos(x 3 , Xj) — 0 — A 31 
cos(x 3 , x 2 ) = 0 = A 32 

and the three-dimensional transformation matrix is 

/ cos 6 sin 6 0 ^ 

A* = — sin d cos 6 0 

V 0 0 1/ 



(1.36b) 



* 3 = *3 




FIGURE 1-10 Coordinate system x lt x 2 , x 3 is rotated an angle 9 ccw about the x 3 -axis. 
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FIGURE 1-1 1 An object undergoes an inversion, which is a reflection about the origin 
of all the axes. 



As a final example, consider the transformation that results in the reflection 
through the origin of all the axes, as in Figure 1-11. Such a transformation is 
called an inversion. In such a case, x[ = — x lt x', = — x 2 , x$ = — x 3 , and 



/- 


-1 


0 


o) 






0 


-1 


0 


(1.37) 


V 


0 


0 


-l) 





In the preceding examples, we defined the transformation matrix A 3 to be 
the result of two successive rotations, each of which was an orthogonal transfor- 
mation: A 3 = A 2 A l . We can prove that the successive application of orthogonal 
transformations always results in an orthogonal transformation. We write 



Xi ijXj, x h 



f^ki x i 



Combining these expressions, we obtain 



X k 5 X j 

= 2 [/* A]fgXj 



Thus, we accomplish the transformation from x t to x" t by operating on x i with the 
{fiA) matrix. The combined transformation will then be shown to be orthogonal 
if (fi A) 1 = (fiA)~ l . The transpose of a product matrix is the product of the 
transposed matrices taken in reverse order (see Problem 1-4); that is, (AB)' = 
B* A*. Therefore 



(/tA)< = Ay 



(1.38) 
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But, because A and fi are orthogonal, A* = A 1 and fj, 1 = /a Multiplying the 
above equation by fiX from the right, we obtain 



( fJL\) t fx\ — X’fJi’fJLX 
= A'lA 
= A 'A 
= 1 

= (M) 'fiX 



Hence 



(mAO = (mA)- 1 (1.39) 

and the fiX matrix is orthogonal. 

The determinants of all the rotation matrices in the preceding examples can 
be calculated according to the standard rule for the evaluation of determinants 
of second or third order: 



[ A | 


= 


An A 12 
A l >i A 2 2 


— Aj ! A 22 


Ai 2 A 2 i 






(1.40) 






An A l2 


^13 












1 A f 


= 


A 2 i A 22 


^23 
















^31 A;n 


A33 














— 


a 22 a 

Ai i 

a 32 a 


23 _ 
33 


a 12 


A 2 i a 23 

^31 ^33 


+ f 

A 3 1 


a 22 

A32 


(1.41) 



where the third-order determinant has been expanded in minors of the first 
row. Therefore, we find, for the rotation matrices used in this section, 

|Aj - |a 2 | - - |a 5 | = l 

but 

I ^6 1 = — 1 

Thus, all those transformations resulting from rotations starting from the original set 
of axes have determinants equal to +1. But an inversion cannot be generated by 
any series of rotations, and the determinant of an inversion matrix is equal to — 1. 
Orthogonal transformations, the determinant of whose matrices is +1, are 
called proper rotations; those with determinants equal to —1 are called im- 
proper rotations. All orthogonal matrices must have a determinant equal to either 
+ 1 or — 1. Here, we confine our attention to the effect of proper rotations and do 
not concern ourselves with the special properties of vectors manifest in improper 
rotations. 
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EXAMPLE 1.4 



Show that | A 2 1 = 1 and | A 6 1 = — 1 . 
Solution. 




1 0 0 

0 0 1 

0 -1 0 




1 

0 



= o — (-1) = 1 





-1 


0 


0 




-1 


0 


1 a 6 I — 


0 


-1 


0 


= -1 


0 


-1 




0 


0 


-1 





= - 1(1 - 0 ) = -1 



1 .8 Definitions of a Scalar and a Vector in Terms 
of Transformation Properties 

Consider a coordinate transformation of the type 



with 



If, under such a transformation, a quantity 4> is unaffected, then <f> is called a 

scalar (or scalar invariant). 

If a set of quantities (A ls A 2 , A 3 ) is transformed from the x, system to the x- 
system by a transformation matrix A with the result 



(1.44) 



_K_ 

II 

-M 

>- 


(1.42) 




(1.43) 




then the quantities A, transform as the coordinates of a point (i.e., according to 
Equation 1.42), and the quantity A = (A lf A 2 , A 3 ) is termed a vector. 



1.9 Elementary Scalar and Vector Operations 

In the following, A and B are vectors (with components A { and B,) and (f > , if/, and 
£ are scalars. 

Addition 



A, + B, = B, + Aj Commutative law (1-45) 

A, + (-B, + C t ) = (A, + B t ) + C, Associative law (1«46) 
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(f> + if/ = (f/ + Commutative law (1*47) 

<p 4- (if/ + £) = (<£ + ijt) + £ Associative law (1 .48) 

Multiplication by a scalar £ 

£A = B is a vector (1.49) 

£<£ = (/' is a scalar (1.50) 

Equation 1 .49 can be proved as follows: 

J J 

= £2a<,a ; = £a; (i.5i) 

and £A transforms as a vector. Similarly, $<p transforms as a scalar. 



1.10 Scalar Product of Two Vectors 



The multiplication of two vectors A and B to form the scalar product is defined 
to be 



AB = Xa ( B, 

i 



(1.52) 



where the dot between A and B denotes scalar multiplication; this operation is 
sometimes called the dot product. 

The vector A has components A u A 2 , A$, and the magnitude (or length) of A 
is given by 

|A| = + Va? + A\ + A§ = A (1.53) 

where the magnitude is indicated by |A| or, if there is no possibility of confu- 
sion, simply by A. Dividing both sides of Equation 1.52 by AB, we have 



A • B ^ A, B , 
AB ~ i A B 



(1.54) 



Aj/A is the cosine of the angle a between the vector A and the xj-axis (see 
Figure 1-12). In general, A,/ A and BJB are the direction cosines A f and A f of 
the vectors A and B: 



= 2a?A f (1.55) 

AB • 

The sum X, A^Af is just the cosine of the angle between A and B (see Equation 

1 . 11 ): 

cos (A, B) = 2A?A? 

* 

or 

r (1.56) 



A • B — AB cos (A, B) 
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*3 




FIGURE 1-12 A vector A is shown in coordinate system x,, x 2 , with its vector 
components A lf A 2 , and A 3 . The vector A is oriented at an angle a 
with the Xj -axis. 



That the product A * B is indeed a scalar may be shown as follows. A and B 
transform as vectors: 



A - - 2 A { jAj, B'j - 2 X ik B k 



(1.57) 



Therefore the product A' * B' becomes 

A'B' = 2 A\B\ 



= 2(2A yA . (2 a,^ 



» \ > 



Rearranging the summations, we can write 

a'-b' = 2(2v*V-«i 

But according to the orthogonality condition, the term in parentheses is just S ]k . 
Thus, 

A'B = 2(2 S ik A:B, 



j \ * 

= 2 AjBj 

= AB 






(1.58) 



Because the value of the product is unaltered by the coordinate transformation, 
the product must be a scalar. 

Notice that the distance from the origin to the point (*], Xq, %) defined by 
the vector A, called the position vector, is given by 



| A | = V / A - A = \/x? + xf + xf — V 2 xf 




CO 



*3 




FIGURE 1-13 The vector A is the position vector of point (x 1; x 2 , x 3 ), and vector B is 
the position vector of point (xj, x 2 , x 3 ) The vector A - B is the 
position vector from (x 1; x 2 , x 3 ) to (x w x 2 , x 3 ). 



Similarly, the distance from the point (xj, x 2 , x 3 ) to another point (xj, x 2 , x 3 ) de- 
fined by the vector B is 

VS(x, - x t f - V(A - B) • (A - B) = |A - B| 

i 

That is, we can define die vector connecting any point with any other point as 
the difference of the position vectors that define the individual points, as in 
Figure 1-13. The distance between the points is then the magnitude of the dif- 
ference vector. And because this magnitude is the square root of a scalar prod- 
uct, it is invariant to a coordinate transformation. This is an important fact and 
can be summarized by the statement that orthogonal transformations are distance- 
preserving transformations. Also, the angle between two vectors is preserved under 
an orthogonal transformation. These two results are essential if we are to suc- 
cessfully apply transformation theory to physical situations. 

The scalar product obeys the commutative and distributive laws: 

A- B = "La i B 1 = S B t A t = B • A (1.59) 

i i 

A- (B + C) = lA t (B + C)i = Ea 2 (A + Q 

i i 

= 2 (A, B i + A, Q - (A-B) + (A-C) (1.60) 



1.11 Unit Vectors 

Sometimes we want to describe a vector in terms of the components along the 
three coordinate axes together with a convenient specification of these axes. For 
this purpose, we introduce unit vectors, which are vectors having a length equal 
to the unit of length used along the particular coordinate axes. For example, the 
unit vector along the radial direction described by the vector R is e R = R/( |R| ). 





24 



1 / MATRICES, VECTORS, AND VECTOR CALCULUS 



There are several variants of the symbols for unit vectors; examples of the most 
common sets are (i, j, k), (e 1? e 2 , e 3 ), (e r , e 9 , e^), and (r, 0, )■ The following 

ways of expressing the vector A are equivalent: 



A = (Aj, A 2 , A 3 ) or A = Ai + e 2 A 2 + e 3 A 3 = 2e,A, 
or A = Aji + A 2 j + A 3 k 



(1.61) 



Although the unit vectors (i, j, k) and (r, 0, <f>) are somewhat easier to use, we 
tend to use unit vectors such as (e^ e 2 , e 3 ) , because of the ease of summation no- 
tation. We obtain the components of the vector A by projection onto the axes: 



A,- = e,-A 



(1.62) 



We have seen (Equation 1.56) that the scalar product of two vectors has a 
magnitude equal to the product of the individual magnitudes multiplied by the 
cosine of the angle between the vectors: 



A • B — AB cos (A, B) 

If any two unit vectors are orthogonal, we have 



(1.63) 



* j 



(1.64) 



Two position vectors are expressed in Cartesian coordinates as A = i -f- 2j — 2k 
and B = 4i + 2j — 3k. Find the magnitude of the vector from point A to point 
B, the angle 6 between A and B, and the component of B in the direction of A. 

Solution. The vector from point A to point B is B — A (see Figure 1-13). 

B - A = 4i + 2j - 3k - (i + 2j - 2k) = 3i - k 
t B — A I — V9 + 1 = VH) 



IB - A | 

From Equation 1.56 



cos 8 = 



A-B (i + 2j - 2k) * (4i + 2j - 3k) 



V9V29 



4 + 4 + 6 

cos 6 = 1 = — = 0.867 

3(V29) 

8 = 30° 



The component of B in the direction of A is B cos 8 and, from Equation 



1.56, 



A-B 14 

B cos 8 = = — = 4.67 

A 3 
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1.12 Vector Product of Two Vectors 

We next consider another method of combining two vectors — the vector prod- 
uct (sometimes called the cross product) . In most respects, the vector product of 
two vectors behaves like a vector, and we shall treat it as such.* The vector prod- 
uct of A and B is denoted by a bold cross X, 

C = A X B (1,65) 

where C is the vector resulting from this operation. The components of C are 
defined by the relation 

(1.66) 

where the symbol e,j k is the permutation symbol or (Levi-Civita density) and has 
the following properties: 

0, if any index is equal to any other index I 

= + 1, if i, j, k form an even permutation of 1, 2, 3 > (1.67) 

— 1, if i, j, k form an odd permutation of 1 , 2, 3 J 

An even permutation has an even number of exchanges of position of two sym- 
bols. Cyclic permutations (for example, 123 — > 231 — > 312) are always even. 
Thus 

e 122 = e 313 = e 211 = 0, etc. 

8 1 23 = e 231 ~ e 312 ~ + 1 
8 132 = e 213 = e 321 = ~ 1 

Using the preceding notation, the components of C can be explicitly evaluated. 
For the first subscript equal to 1, the only nonvanishing s ijk are e 123 and e 132 — 
that is, for j, k = 2, 3 in either order. Therefore 

Q = 2>£j;*A :B k = £j 2 gA 2 .B 3 + £ J32 A 3 fi 2 

],k 1 J 

= A 2 B 3 A 3 B 2 (1.68a) 

Similarly, 

C 2 = A 3 Bj - AjS 3 (1.68b) 

C 3 = A l B 2 - A 2 B, (1.68c) 

Consider now the expansion of the quantity [ABsin(A, B)] 2 = (ABsinfl) 2 : 

A 2 B 2 sin 2 0 = A 2 B 2 - A 2 B 2 cos 2 0 

= ( A 2 B s - A 3 £ 2 ) 2 + (A 3j B, - A,B 3 ) 2 + {A x B 2 - A 2 B x ) 2 (1.69) 




*The product actually produces an axial vector, but the term vector product is used to be consistent 
with popular usage. 
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A 

FIGURE 1-14 The magnitude of the vector C determined by C = A X B has a 

magnitude given by the area of the parallelogram AB sin d, where 6 
is the angle between the vectors A and B. 

where the last equality requires some algebra. Identifying the components of C 
in the last expression, we can write 

(AB sin 0) 2 = C? + C§ + C| = |C 2 | = C 2 (1.70) 

If we take the positive square root of both sides of this equation, 

C = AB sin 8 (1.71) 

This equation states that if C = A x B, the magnitude of C is equal to the prod- 
uct of the magnitudes of A and B multiplied by the sine of the angle between 
them. Geometrically, AB sin 6 is the area of the parallelogram defined by the 
vectors A and B and the angle between them, as in Figure 1-14. 



EXAMPLE 1.6 



Show by using Equations 1.52 and 1.66 that 

A- (B X D) = D • (A x B) (1.72) 

Solution. Using Equation 1.66, we have 

(B x D), — 2 e,j k BjD k 
Using Equation 1.52, we have 

A* (B x D) = '2e l)k A l B j D k (1.73) 

Similarly, for the right-hand side of Equation 1.72, we have 

D • (A x B) = 2e*A AjB k 

i,j,k j 

From the definition (Equation 1.67) of e yi , we can interchange two adjacent in- 
dices of e y *, which changes the sign. 

D • (A X B) = 2 — e )lk D t A B k 

l t J,k 

~ 2* )kt A A B t (1.74) 

i ,j,k J J 

Because the indices i, j, k are dummy and can be renamed, the right-hand sides 
of Equations 1.73 and 1.74 are identical, and Equation 1.72 is proved. Equation 
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1.72 can also be written as A ■ (B X D) = (A X B) ■ D, indicating that the scalar 
and vector products can be interchanged as long as the vectors stay in the order 
A, B, D. Notice that, if we let B — A, we have 

A • (A X D) = D • (A X A) = 0 

showing that A X D must be perpendicular to A. 



A X B (i.e., C) is perpendicular to the plane defined by A and B because 
A • (A X B) = 0 and B • (A X B) = 0. Because a plane area can be represented 
by a vector normal to the plane and of magnitude equal to the area, C is evi- 
dently such a vector. The positive direction of C is chosen to be the direction of 
advance of a right-hand screw when rotated from A to B. 

The definition of the vector product is now complete; components, magni- 
tude, and geometrical interpretation have been given. We may therefore reason- 
ably expect that C is indeed a vector. The ultimate test, however, is to examine 
the transformation properties of C, and C does, in fact, transform as a vector 
under a proper rotation. 

We should note the following properties of the vector product that result 



from the definitions: 

(a) 


A X B = -B X A 


(1.75) 


but, in general, 

(b) 


Ax (B x C) / (A x B) x C 


(1.76) 


Another important result (see Problem 1-22) is 

A x (B x C) = (A • C)B - (A • B)C 


(1.77) 



EXAMPLE 1.7 



Find the product of (A X B) ■ (C X D). 

Solution. 

(A x B), - e ljk AjB k 

),k j 

(CXD),. = S e ilm C l D m 

The scalar product is then computed according to Equation 1 .52: 

Rearranging the summations, we have 

(AXB).(CXD)^ (IS A ] B k C,D m 

where the indices of the e’s have been permuted (twice each so that no sign 
change occurs) to place in the third position the index over which the sum is 
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carried out. We can now use an important property of the e i;k (see Problem 1-22): 

(1.78) 

We therefore have 

(A X B) • (C X D) =2 - S Jm S u )A J B t C,D m 

t,m 

Carrying out the summations over j and k, the Kronecker deltas reduce the ex- 
pression to 

(A X B) ■ (C X D) " 2 (A l B m C l D m - A m B,C l D m ) 

This equation can be rearranged to obtain 

(A X B) • (C X D) = feA,C,Y2B«Bj - (£b,C,Y2a„.D,) 




Because each term in parentheses on the right-hand side is just a scalar product, 
we have, finally, 

(A x B) • (C X D) = (A-C)(B-D) - (B-C)(A-D) 



The orthogonality of the unit vectors e* requires the vector product to be 



e { x = e k i,j, kin cyclic order (1.79a) 

We can now use the permutation symbol to express this result as 

e< X (1.79b) 

The vector product C = A x B, for example, can now be expressed as 

C = ^e ijk e,A j B k (1.80a) 

l,j,k J 

By direct expansion and comparison with Equation 1.80a, we can verify a de- 
terminantal expression for the vector product: 

ei e 2 €3 

C = A X B = A, A 2 A 3 (1.80b) 

We state the following identities without proof: 

A • (B X C) — B • (C X A) = C • (A x B) = ABC (1.81) 

A x (B x C) = (A*C)B — (A * B)C (1.82) 



(A x B) • (C x D) = A- [B x (C x D)] 

= A- [(B-D)C - (B • C)D] 

= (A* C) (B • D) - (A*D) (B * C) 

(A x B) x (C x D) = [(A x B) • D]C - [(A x B) • C]D 

= (ABD)C - (ABC)D = (ACD)B - (BCD)A 



(1.83) 

(1.84) 
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1.13 Differentiation of a Vector 
with Respect to a Scalar 



If a scalar function cf) = </>(s) is differentiated with respect to the scalar variable s, 
then, because neither part of the derivative can change under a coordinate 
transformation, the derivative itself cannot change and must therefore be a 
scalar; that is, in the x, and x' coordinate systems, <f> = <j>' and s = s', so d<f> = d<p' 
and ds = ds'. Hence 

d<f> d(f>' ( d<p\ 

ds ds' \ds J 



Similarly, we can formally define the differentiation of a vector A with re- 
spect to a scalar s. The components of A transform according to 

A’ = 2 \ tj A, (1.85) 

Therefore, on differentiation, we obtain (because the A y are independent of s') 



dA’j 

ds’ 



ds' ^ XijAj ^ A * 



dA ! 

ds' 



Because s and s' are identical, we have 

dA \ _ ( dA \ dA j \ 

ds' \ds) .7 ij \ ds ) 

Thus the quantities dAj/ds transform as do the components of a vector and 
hence are the components of a vector, which we can write as dX/ ds. 

We can give a geometrical interpretation to the vector dX/ ds as follows. First, 
for dX/ ds to exist, A must be a continuous function of the variable s: A = A(s). 
Suppose this function is represented by the continuous curve T in Figure 1-15; at 
the point P, the variable has the value s, and at Qit has the value s + As. The de- 
rivative of A with respect to s is then given in standard fashion by 

^ = lim - lim (1.86a) 

ds As— *•() As As— >0 As 



*2 




FIGURE 1-15 The vector A(s) traces out the function T(s) as the variable s changes. 
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The derivatives of vector sums and products obey the rules of ordinary vec- 
tor calculus. For example, 



d 

-{A + B) 
ds 



d 

ds 



(A-B) 



_rf 

ds 



(A x B) 




dA dB 

ds ds 


(1.86b) 


dB dA 
— A • 1 • B 

ds ds 


(1.86c) 


dB dA 

= A x — + — x B 

ds ds 


(1.86d) 


, dA dxf) 

— 4> 1 A 

ds ds 


(1.86e) 



and similarly for total differentials and for partial derivatives. 



1.14 Examples of Derivatives — 

Velocity and Acceleration 

Of particular importance in the development of the dynamics of point particles 
(and of systems of particles) is the representation of the motion of these parti- 
cles by vectors. For such an approach, we require vectors to represent the posi- 
tion, velocity, and acceleration of a given particle. It is customary to specify the 
position of a particle with respect to a certain reference frame by a vector r, which 
is in general a function of time: r — r(f). The velocity vector v and the acceleration 
vector a are defined according to 



v = 



dr 

dt 



— r 



dv d 2 r 

dt dt 2 



( 1 . 87 ) 

( 1 . 88 ) 



where a single dot above a symbol denotes the first time derivative, and two dots 
denote the second time derivative. In rectangular coordinates, the expressions 
for r, v, and a are 



r = Xjej 4- x 2 e 2 -t- x 3 e 3 = 2 x,e, 

l 

V . v dxi 
v = r = ZjX:e, — Zj — 

• ‘dt 



— v = r = e, = 2 



d 2 x, 
dt 2 



e, 



Position 

Velocity 

Acceleration 



( 1 . 89 ) 



Calculating these quantities in rectangular coordinates is straightforward because 
the unit vectors are constant in time. In nonrectangular coordinate systems, 
however, the unit vectors at the position of the particle as it moves in space are 
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not necessarily constant in time, and the components of the time derivatives of r 
are no longer simple relations, as in Equation 1.89. We do not discuss general 
curvilinear coordinate systems here, but plane polar coordinates, spherical coordi- 
nates, and cylindrical coordinates are of sufficient importance to warrant a discus- 
sion of velocity and acceleration in these coordinate systems.* 

To express v and a in plane polar coordinates, consider the situation in 
Figure 1-16. A point moves along the curve s(t) and in the time interval 
t 2 ~ t x = dt moves from P (l) to P (2) . The unit vectors, e r and e 0 , which are or- 
thogonal, change from e£ 1} to e) 2) and from e^ 1 ' to e< 2> . The change in e, is 



e< 2) — ej- 15 — de r 


(1.90) 


which is a vector normal to e T (and, therefore, in 
the change in e e is 


the direction of e e ) . Similarly, 


e (2) - e<h = de 6 


(1.91) 


which is a vector normal to e e . We can then write 




de r = dOe e 


(1.92) 


and 




de 9 = —dde r 


(1.93) 



where the minus sign enters the second relation because de e is directed opposite 
to e r (see Figure 1-16). 



FIGURE 1-16 



s(t) 




An object traces out the curve $(t) over time. The unit vectors e r and e 0 
and their differentials are shown for two position vectors rj and r 2 . 



* Refer to the figures in Appendix F for the geometry of these coordinate systems. 
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Equations 1.92 and 1.93 are perhaps easier to see by referring to Figure 1-16. 
In this case, de T subtends an angle dd with unit sides, so it has a magnitude of dd. 
It also points in the direction of e 6 , so we have de T = dde 9 . Similarly, de 6 subtends 
an angle dd with unit sides, so it also has a magnitude of dd, but from Figure 1-16 
we see that de g points in the direction of — e r so we have de e = — dde r . 

Dividing each side of Equations 1.92 and 1.93 by dt, we have 

(1.94) 

(1.95) 

If we express v as 



e r = de e 
= ~de r 



v 



dr 

dt 



d L 

dt 



(re r ) 



fe r + re r 



(1.96) 



we have immediately, using Equation 1 .94, 



v = r = fe r + rde e 



(1.97) 



so that the velocity is resolved into a radial component r and an angular (or 
transverse) component rd. 

A second differentiation yields the acceleration: 



a = ^(*e r + r ® e o) 

= re r + fe r + rde 9 + rde e + rde e 
= (r — rd 2 )e r + (rd + 2rd)e e > (1.98) 



so that the acceleration is resolved into a radial component (r — rd 2 ) and an 

4» • 

angular (or transverse) component (rd + 2 rd). 

The expressions for ds, ds 2 , v 2 , and v in the three most important coordi- 
nate systems (see also Appendix F) are 



Rectangular coordinates (x, y, z) 

ds = dxje i + dx 2 e 2 + dx 3 eg| 
ds 2 = dx i + dx \ + dx 1 
v 2 = X? + x| + x§ 
v = xjej + x 2 e 2 + x 3 e 3 

Spherical coordinates (r, d, </>) 

ds = dre r + rdde e + r sin d d(f>e$ | 
ds 2 = dr 2 + r 2 dd 2 + r 2 sin 2 0 d<f> 2 
v 2 — r 2 + r 2 d 2 + r 2 sin 2 0 <j> 2 
v = fe r + rde e + r sin d 






( 1 . 100 ) 



(1.99) 
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{The expressions for plane polar coordinates result from Equation 1.100 by set- 
ting dcf> - 0.) 



Cylindrical coordinates (r, d>, z) 

ds — dre T + rdcfre# + dze, 
ds 2 = dr 2 + r 2 d(f> 2 + dz? 



v 



2 — r 2 + r 2 <j> 2 + z 2 



v = fe r T + ze z 



( 1 . 101 ) 



EXAMPLE 1.8 



Find the components of the acceleration vector a in cylindrical coordinates. 

Solution. The velocity components in cylindrical coordinates were given in 
Equation 1.101. The acceleration is determined by taking the time derivative of v. 

d d,, 

a = — v = — (re r + + ze z ) 

dt dt 

- re r + re r + + ze z + ze z 

We need to find the time derivative of the unit vectors e,, e^,, and e 2 . The 
cylindrical coordinate system is shown in Figure 1-17, and in terms .of the (x, y, z) 
components, the unit vectors e r , e^, and e z are 

e r = ( cos <f>, sin <f>, 0) 
e $ — ( — sin (j>, cos <£, 0) 
e, = (0, 0, 1) 



z 




FIGURE 1-17 The cylindrical coordinate system (r, < f>, z) are shown with respect to the 
Cartesian system (x, y, z). 
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The time derivatives of the unit vectors are found by taking the derivatives of 
the components. 

e r = (— <j> sin 4>, <p cos <p, 0) = — 

= ( — (f> cos sin (f>, 0) = — (j>e r 

e, - 0 

We substitute the unit vector time derivatives into the above expression for a. 

• • M • 

a = re r + rtpe^ + rcfre# + r<£e^ — r(f> z e r + ze 2 
= (r — r<£ 2 )e r + {r<j> + 2r4>)e tf) + ze. 



1.15 Angular Velocity 

A point or a particle moving arbitrarily in space may always be considered, at a 
given instant , to be moving in a plane, circular path about a certain axis; that is, 
the path a particle describes during an infinitesimal time interval 8t may be rep- 
resented as an infinitesimal arc of a circle. The line passing through the center 
of the circle and perpendicular to the instantaneous direction of motion is 
called the instantaneous axis of rotation. As the particle moves in the circular 
path, the rate of change of the angular position is called the angular velocity: 

d0 . 

a) = ~ = 6 (1.102) 

at 

Consider a particle that moves instantaneously in a circle of radius R about 
an axis perpendicular to the plane of motion, as in Figure 1-18. Let the position 
vector r of the particle be drawn from an origin located at an arbitrary point O 
on the axis of rotation. The time rate of change of the position vector is the 
linear velocity vector of the particle, r — v. For motion in a circle of radius R, the 
instantaneous magnitude of the linear velocity is given by 

dd 

v - R— = R(o (1.103) 

(Xt 

The direction of the linear velocity v is perpendicular to r and in the plane of the 
circle. 

It would be very convenient if we could devise a vector representation of 
the angular velocity (say, to) so that all the quantities of interest in the motion 
of the particle could be described on a common basis. We can define a direction 
for the angular velocity in the following manner. If the particle moves instanta- 
neously in a plane, the normal to that plane defines a precise direction in 
space— or, rather — two directions. We may choose as positive that direction correspon- 
ding to the direction of advance of a right-hand screw when turned in the same 
sense as the rotation of the particle (see Figure 1-18). We can also write the mag- 
nitude of the linear velocity by noting that R = rsina. Thus 



v — ra> sin a 



(1.104) 
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FIGURE 1-18 A particle moving ccw about an axis according to the right-hand rule 
has an angular velocity to = v X r about that axis. 



Having defined a direction and a magnitude for the angular velocity, we note 
that if we write 



v = to X r 



(1.105) 



then both of these definitions are satisfied, and we have the desired vector rep- 
resentation of the angular velocity. 

We should note at this point an important distinction between finite and in- 
finitesimal rotations. An infinitesimal rotation can be represented by a vector 
(actually, an axial vector) , but a finite rotation cannot. The impossibility of de- 
scribing a finite rotation by a vector results from the fact that such rotations do 
not commute (see the example of Figure 1-9), and therefore, in general, differ- 
ent results will be obtained depending on the order in which the rotations are 
made. To illustrate this statement, consider the successive application of two fi- 
nite rotations described by the rotation matrices A, and A 2 . Let us associate the 
vectors A and B in a one-to-one manner with these rotations. The vector sum is C = 
A -f B, which is equivalent to the matrix A 3 = A 2 A ] . But because vector addition 
is commutative, we also have C = B + A, with A 4 = AjA 2 . But we know that 
matrix operations are not commutative, so that in general A 3 =£ A 4 . Hence, the 
vector C is not unique, and therefore we cannot associate a vector with a finite 
rotation. 

Infinitesimal rotations do not suffer from this defect of noncommutation. We 
are therefore led to expect that an infinitesimal rotation can be represented by a 
vector. Although this expectation is, in fact, fulfilled, the ultimate test of the vec- 
tor nature of a quantity is contained in its transformation properties. We give 
only a qualitative argument here. 

Refer to Figure 1-19. If the position vector of a point changes from r to r + 
8r, the geometrical situation is correctly represented if we write 



8r = 50 x r 



(1.106) 
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50 




FIGURE 1-19 The position vector r changes to r + 5r by an infinitesimal 
rotation angle 80. 

f 

where 50 is a quantity whose magnitude is equal to the infinitesimal rotation f 

angle and that has a direction along the instantaneous axis of rotation. The j 

mere fact that Equation 1.106 correctly describes the situation illustrated in I 

Figure 1-19 is not sufficient to establish that 60 is a vector. (We reiterate that the • j 

true test must be based on the transformation properties of 60.) But if we show j 

that two infinitesimal rotation “vectors” — 60! and 602 — actually commute, the sole « 

objection to representing a finite rotation by a vector will have been removed. f 

Let us consider that a rotation 60, takes r into r + 6r,, where 6r, = 60 t X r. f 

If this is followed by a second rotation 60 2 around a different axis, the initial po- ■ 

sition vector for this rotation i$ r 4- 5r,. Thus \ 

6r 2 = S0 2 x (r + 6r,) 

and the final position vector for 60! followed by 50 2 is j 

r + 5r 12 — r + [50, X r + 60 2 X (r + 5r,)] 

Neglecting second-order infinitesimals, then, 

6r 12 = 60, x r + 60 2 x r (1.107) 

Similarly, if 60 2 is followed by 60,, we have 

r + 6r 21 = r + [60 2 x r + 50, x (r + 5r 2 )] 
or 

6r 21 = 50 2 x r + 50, x r (1.108) 

Rotation vectors 6r, 2 and 5r 2 , are equal, so the rotation “vectors” 50, and 50 2 do 
commute. It therefore seems reasonable that 60 in Equation 1.106 is indeed a 
vector. 
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It is the fact that 50 is a vector that allows angular velocity to be represented 
by a vector, because angular velocity is the ratio of an infinitesimal rotation angle 
to an infinitesimal time: 



50 

to = — • 

8t 

Therefore, dividing Equation 1.106 by 8t, we have 

Sr 60 

— — — x r 
8t St 

or, in passing to the limit, 8t — > 0, 

v — to x r 



as before. 



1.16 Gradient Operator 

We now turn to the most important member of a class called vector differential 
operators — the gradient operator. 

Consider a scalar <f> that is an explicit function of the coordinates x ; and, 
moreover, is a continuous, single-valued function of these coordinates through- 
out a certain region of space. Under a coordinate transformation that carries the 
x, into the x', 4)'(x[, x' 2 , x 3 ) = x 2 , x 3 ), and by the chain rule of differentia- 

tion, we can write 



()<f) ' ^ d<b dXj 

dx[ j dXjdx[ 

The case is similar for d4>'/dx 2 and d<t>'/dx 3 , so in general we have 

d(f> ' ^ d(j> dx. 



uyj '"y 

dx’i j dxj dx[ 

The inverse coordinate transformation is 

Xj 

Differentiating, 



dxj 



dxi dx;(? A **‘, 






But the term in the last parentheses is just 5 rt , so 



dxj 



(1.109) 



( 1 . 110 ) 



( 1 . 111 ) 



( 1 . 112 ) 
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Substituting Equation 1.113 into Equation 1.110, we obtain 



— = Xa — 

ax; j ij d Xj 



(1.114) 



Because it follows the correct transformation equation of a vector (Equation 
1.44), the function d0/dx ; is the jth component of a vector termed the gradient 
of the function 0. Note that even though 0 is a scalar, the gradient of 0 is a vector. 
The gradient of 0 is written either as grad 0 or as V0(“del” 0). 

Because the function 0 is an arbitrary scalar function, it is convenient to de- 
fine the differential operator described in the preceding in terms of the gradient 
operator : 



(grad), = V, = 



d 

dx, 



(1.115) 



We can express the complete vector gradient operator as 

d 

grad = V = Ze, — | Gradient 

1 ox. 



(1.116) 



The gradient operator can (a) operate directly on a scalar function, as in 
V0; (b) be used in a scalar product with a vector function, as in V ■ A (the diver- 
gence (div) of A); or (c) be used in a vector product with a vector function, as in 
V x A (the curl of A). We present the grad, divergence, and curl: 



V' d0 

grad0 = V0 = ^e, — (1.117a) 

i ()X, 

•sr' dAj 

div A = V * A = aL — - (1.117b) 

i iix l 

V dA* 

curl A — V X A = e ! (1.117c) 

i,j,k J dXj 

To see a physical interpretation of the gradient of a scalar function, consider 
the three-dimensional and topographical maps of Figure 1-20. The closed loops 
of part b represent lines of constant height. Let 0 denote the height at any point 
0 = 0(xj, x 3 ). Then 

d0 = 2 — dx, = 2(V0),dx, 

i OX, i 



The components of the displacement vector ds are the incremental displace- 
ments in the direction of the three orthogonal axes: 

ds — (dX], dx 2 , dx 3 ) (1.118) 



Therefore 



d0 = (V0) • ds 



(1.119) 
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y 




(b) 

FIGURE 1-20 (a) A three-dimensional contour map can be represented by (b) a 

topographical map of lines <f> representing constant height. The 
gradient V(f> represents the direction perpendicular to the constant 
<f> lines. 



Let ds be directed tangentially along one of the isolatitude lines (i.e., along 
a line for which </> = const.), as indicated in Figure 1-20. Because 0 — const, for 
this case, d(f) = 0. But, because neither V 4> nor ds is in general zero, they must there- 
fore be perpendicular to each other. Thus V<£ is normal to the line (or in three 
dimensions, to the surface) for which </> = const. 
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The maximum value of d4> results when V0 and ds are in the same direc- 
tion; then, 

(#)max = |V0|di, for V<£||ds 



or 




Therefore, Vcf) is in the direction of the greatest change in <fi. 
We can summarize these results as follows: 



( 1 . 120 ) 



1. The vector is, at any point, normal to the lines or surfaces for which 0 = 
const. 

2. The vector V<£ has the direction of the maximum change in 4>. 

3. Because any direction in space can be specified in terms of the unit vector n 
in that direction, the rate of change of < p in the direction of n (the directional 
derivative of (f>) can be found from n ■ V(fi = d(f>/dn. 



The successive operation of the gradient operator produces 

_ _ XT' d d v tl 2 

V • V = 2j- r — —z 

i dXj dx, ‘ ox i 



( 1 . 121 ) 



This important product operator, called the Laplacian,* is also written 




When the Laplacian operates on a scalar, we have, for example, 



Vfy = 



"S — 

i 



( 1 . 122 ) 



( 1 . 123 ) 



1.17 Integration of Vectors 



The vector resulting from the volume integration of a vector function A = A(x,) 
throughout a volume V is given byf 



A dv 
. v 




A 2 dv, 
. v 




( 1 . 124 ) 



♦After Pierre Simon Laplace (1749-1827); the notation V 2 is ascribed to Sir William Rowan 
Hamilton. 

’The symbol J v actually represents a triple integral over a certain volume V. Similarly, the symbol f s 
stands for a double integral over a certain surface 5. 
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FIGURE 1-21 The differential da is an element of area of the surface. Its direction is 
normal to the surface. 



Thus, we integrate the vector A throughout V simply by performing three sepa- 
rate, ordinary integrations. 

The integral over a surface S of the projection of a vector function A = A(x ; ) 
onto the normal to that surface is defined to be 



J A* da 

where da is an element of area of the surface (Figure 1-21). We write d a as a vec- 
tor quantity because we may attribute to it not only a magnitude da but also a di- 
rection corresponding to the normal to the surface at the point in question. If 
the unit normal vector is n, then 

da — n da (1.125) 

Thus, the components of da are the projections of the element of area on the 
three mutually perpendicular planes defined by the rectangular axes: 



Therefore, we have 



or 



= dx 2 dx$, etc. 



j A - da 



J A* n da 



A* da — 



2a da. 



Js 



(1,126) 



(1.127) 



(1.128) 



Equation 1.127 states that the integral of A over the surface S is the integral of 
the normal component of A over this surface. 

The normal to a surface may be taken to lie in either of two possible direc- 
tions (“up” or “down”); thus the sign of n is ambiguous. If the surface is dosed , we 
adopt the convention that the outward normal is positive. 

The line integral of a vector function A = A(Xj) along a given path extend- 
ing from the point B to the point C is given by the integral of the component of 
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C 




FIGURE 1-22 The element ds is an element of length along the given path from B to 
C. Its direction is along the path at a given point. 

A along the path 



1 A - ds 

JHC 



2 A, dx, 



JHC 1 



(1.129) 



The quantity ds is an element of length along the given path (Figure 1-22). The 
direction of ds is taken to be positive along the direction the path is traversed. In 
Figure 1-22 at point P, the angle between ds and A is less than 7r/2, so A • ds is 
positive at this point. At point the angle is greater than 7 t/ 2, and the contri- 
bution to the integral at this point is negative. 

It is often useful to relate certain surface integrals to either volume integrals 
(Gauss’s theorem) or line integrals (Stokes’s theorem). Consider Figure 1-23, 
which shows a closed volume ^enclosed by the surface S. I.et the vector A and its 
first derivatives be continuous throughout the volume. Gauss’s theorem states 
that the surface integral of A over the closed surface S is equal to the volume in- 
tegral of the divergence of A (V • A) throughout the volume V enclosed by the 
surface S. We write this mathematically as 



J A • da. 



J V • A dv 



(1.130) 



Gauss’s theorem is sometimes also called the divergence theorem. The theorem is 
particularly useful in dealing with the mechanics of continuous media. 

See Figure 1-24 for the physical description needed for Stokes’s theorem, 
which applies to an open surface S and the contour path C that defines the sur- 
face. The curl of the vector A (V x A) must exist and be integrable over the en- 
tire surface S. Stokes’s theorem states that the line integral of the vector A 
around the contour path C is equal to the surface integral of the curl of A over 
the surface defined by C. We write it mathematically as 



r 



A* ds 

. c 



(V x A) • da. 

.s 



(1.131) 



where the line integral is around the closed contour path C. Stokes’s theorem is 
particularly useful in reducing certain surface integrals (two dimensional) to, it 
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Surface S 



FIGURE 1-23 The differential da is an element of area on a surface S that surrounds 
a closed volume V. 




Surface S 



Contour C 



FIGURE 1-24 A contour path C defines an open surface S. A line integral around the 
path C and a surface integral over the surface 5 is required for Stokes’s 
theorem. 



is hoped, a simpler line integral (one dimensional). Both Gauss’s and Stokes’s 
theorems have wide application in vector calculus. In addition to mechanics, 
they are also useful in electromagnetic applications and in potential theory. 



PROBLEMS 



1-1. Find the transformation matrix that rotates the axis jc a of a rectangular coordinate 
system 45 ° toward x-i around the A^-axis. 

1-2. Prove Equations 1.10 and 1.11 from trigonometric considerations. 

1-3. Find the transformation matrix that rotates a rectangular coordinate system 
through an angle of 120° about an axis making equal angles with the original three 
coordinate axes. 

1-4. Show 

(a) (AB)' = BA' (b) (AB ) -1 = B 'A ' 

1-5. Show by direct expansion that |A| 2 = 1. For simplicity, take A to be a two- 
dimensional orthogonal transformation matrix. 



44 



1 / MATRICES, VECTORS, AND VECTOR CALCULUS 



1-6. Show that Equation 1.15 can be obtained by using the requirement that the trans- 
formation leaves unchanged the length of a line segment. 

1-7. Consider a unit cube with one corner at the origin and three adjacent sides lying 
along the three axes of a rectangular coordinate system. Find the vectors describ- 
ing the diagonals of the cube. What is the angle between any pair of diagonals? 

1-8. Let A be a vector from the origin to a point P fixed in space. Let r be a vector from 
the origin to a variable point (?(x 1; x 2 , x 3 ) . Show that 

A • r — A 2 

is the equation of a plane perpendicular to A and passing through the point P. 

1-9. For the two vectors 

A = i + 2j - k, B = -2i + 3j + k 

find 

(a) A — B and |A — B| (b) component of B along A (c) angle between A and B 
(d)AXB (e) (A - B) X (A + B) 

1-10. A particle moves in a plane elliptical orbit described by the position vector 

r = 2b sin (oti + b cos cot j 

(a) Find v, a, and the particle speed. 

(b) What is the angle between v and a at time l = 7t/2o)? 

1-11. Show that the triple scalar product (A X B) • C can be written as 

Aj A 2 A 3 
(A x B) • C = B x A, 

Q Q 

Show also that the product is unaffected by an interchange of the scalar and vector 
product operadons or by a change in the order of A, B, C, as long as they are in 
cyclic order; that is, 

(A x B) • C = A- (B x C) = B • (C x A) = (C x A) • B, etc. 

We may therefore use the notation ABC to denote the triple scalar product. Finally, 
give a geometric interpretation of ABC by computing the volume of the paral- 
lelepiped defined by the three vectors A, B, C. 

1-12. Let a, b, c be three constant vectors drawn from the origin to the points A, B, C. 
What is the distance from the origin to the plane defined by the points A, B, C? 
What is the area of the triangle ABC ? 

1-13. X is an unknown vector satisfying the following relations involving the known vec- 
tors A and B and the scalar 4>, 

A X X = B, A X = <£. 

Express X in terms of A, B, <f>, and the magnitude of A. 
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1-14. Consider the following matrices: 

(l 2 -1\ / 2 1 0\ / 2 1\ 

A = 0 3 1, B- 0 -1 2, C = 4 3 

\2 01/ \ 1 1 3/ \l 0/ 

Find the following 

(a) | AB | (b) AC (c) ABC (d) AB - B'A' 

1-15. Find the values of a needed to make the following transformation orthogonal. 

(' 0 °\ 

, I 0 a -a I 

\0 a a/ 

1-16. What surface is represented by r • a = const, that is described if a is a vector of con- 
stant magnitude and direction from the origin and r is the position vector to the 
point P(x i, x 2 , *s) on surface? 

1-17. Obtain the cosine law of plane trigonometry by interpreting the product (A — B) • 
(A - B) and the expansion of the product. 

1-18. Obtain the sine law of plane trigonometry by interpreting the product A X B and 
the alternate representation (A — B) X B. 

1-19. Derive the following expressions by using vector algebra: 

(a) cos {a — 1 8 ) — cos a cos /3 + sin or sin /3 

(b) sin (a ~ fi) — sin a cos (3 — cos a sin fi 

1-20. Show that 

(a) 2e (> * 8y = 0 (b) e tjk = 2 8 U (c) e ijk = 6 

1-21. Show (see also Problem 1-11) that 

ABC = 2 A; B ] C k 

1-22. Evaluate the sum (which contains 3 terms) by considering the result for 

all possible combinations of i, j, l, m\ that is, 

(a) i ~ j (b) i = l (c ) i — m (d) j = l (e) j= m (f ) 1= m 

(g) i ¥= l or m (h) j # lor m 

| Show that 

; ^ E !,k E lmk ~ ~ Sim&jl 

s. * 

I) 

\ and then use this result to prove 



A X (B X C) — (A * C)B — (A • B)C 
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1-23. Use the s tjk notation and derive the identity 

(A X B) X (C x D) = (ABD)C - (ABC)D 

1-24. Let A be an arbitrary vector, and let e be a unit vector in some fixed direction. Show 
that 



A = e(A ■ e) + e X (A X e) 

What is the geometrical significance of each of the two terms of the expansion? 

1-25. Find the components of the acceleration vector a in spherical coordinates. 

1-26. A particle moves with v = const, along the curve r—k{ 1 + cos d) (a cardioid). Find 
r • e, = a • e r , | a | , and 9. 

1-27. If r and f = v are both explicit functions of time, show that 

— [r X (v X r)] = r 2 a + (r-v)v — (v 2 + r-a)r 
dt 



1-28. Show that 



V(ln | r | ) = — 
r £ 

1-29. Find the angle between the surfaces defined by r 2 = 9 and x + y + z 2 — 1 at the 
point (2, —2, 1). 

1-30. Show that V(<j> ip) = 4>Vtp + 

1-31. Show that 



(a) Vr" = nr (M 2) r (b) V/(r) 



rdf 

r dr 



(c) V 2 (ln r) — 



1 



1-32. Show that 




ar-T + 2 br • r )dt — ar 2 + br 2 + const. 



where r is the vector from the origin to the point (*], x 2 , x 3 ). The quantities rand r 
are the magnitudes of the vectors r and r, respectively, and a and b are constants. 

1-33. Show that 




where C is a constant vector. 
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1-34. Evaluate the integral 

A x A dt 

j 

1-35. Show that the volume common to the intersecting cylinders defined by x 2 4 - y 2 = a 2 
and x 2 + z 2 = a 2 is V = 16a 3 /3. 

1-36. Find the value of the integral A ■ da, where A = xi — yj + zk and S is the closed 
surface defined by the cylinder c 2 = x 2 + y 2 . The top and bottom of the cylinder 
are at z = d and 0, respectively. 

1-37. Find the value of the integral / S A> da, where A — (x 2 + y 2 + z 2 ) (xi + yj + zk) and 
the surface S is defined by the sphere R 2 = x 2 + y 2 + z 2 . Do the integral directly 
and also by using Gauss’s theorem. 

1-38. Find the value of the integral Js (V X A) • da if the vector A = yi + zj + xk and 5 is 
the surface defined by the paraboloid z = 1 — x 2 — y 2 , where z > 0. 

1-39. A plane passes through the three points (x, y, z) = (1, 0, 0), (0, 2, 0), (0, 0, 3). 

(a) Find a unit vector perpendicular to the plane, (b) Find the distance from the 
point (1, 1, 1) to the closest point of the plane and the coordinates of the closest 
point. 

1-40. The height of a hill in meters is given by z = 2 xy — 3X 2 — 4 f — 18x + 28y + 12, 
where x is the distance east and y is the distance north of the origin, (a) Where is 
the top of the hill and how high is it? (b) How steep is the hill at x = y = 1, that is, 
what is the angle between a vector perpendicular to the hill and the z axis? (c) In 
which compass direction is the slope at x = y = 1 steepest? 

1-41. For what values of a are the vectors A = 2oi — 2j + <zk and B = ai + 2aj + 2k 
perpendicular? 




CHAPTER 

Newtonian Mechanics — 

Single Particle 



2.1 Introduction 

The science of mechanics seeks to provide a precise and consistent descrip- 
tion of the dynamics of particles and systems of particles, that is, a set of phys- 
ical laws mathematically describing the motions of bodies and aggregates of 
bodies. For this, we need certain fundamental concepts such as distance and 
time. The combination of the concepts of distance and time allows us to 
define the velocity and acceleration of a particle. The third fundamental 
concept, mass, requires some elaboration, which we give when we discuss 
Newton’s laws. 

Physical laws must be based on experimental fact. We cannot expect a pri- 
ori that the gravitational attraction between two bodies must vary exactly as 
the inverse square of the distance between them. But experiment indicates 
that this is so. Once a set of experimental data has been correlated and a pos- 
tulate has been formulated regarding the phenomena to which the data refer, 
then various implications can be worked out. If these implications are all veri- 
fied by experiment, we may believe that the postulate is generally true. The 
postulate then assumes the status of a physical law. If some experiments dis- 
agree with the predictions of the law, the theory must be modified to be con- 
sistent with the facts. 

Newton provided us with the fundamental laws of mechanics. We state these 
laws here in modem terms, discuss their meaning, and then derive the implications 
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of the laws in various situations.* But the logical structure of the science of me- 
chanics is not straightforward. Our line of reasoning in interpreting Newton’s laws 
is not the only one possible.* We do not pursue in any detail the philosophy of me- 
chanics but rather give only sufficient elaboration of Newton’s laws to allow us to 
continue with the discussion of classical dynamics. We devote our attention in this 
chapter to the motion of a single particle, leaving systems of particles to be dis- 
cussed in Chapters 9 and 11-13. 



2.2 Newton’s Laws 

We begin by simply stating in conventional form Newton’s laws of mechanics*: 

I. A body remains at rest or in uniform, motion unless acted upon by a force. 

II. A body acted upon by a force moves in such a manner r that the time rate of change of 
momentum equals the force. 

III. If two bodies exert forces on each other, these forces are equal in magnitude and oppo- 
site in direction. 

These laws are so familiar that we sometimes tend to lose sight of their true 
significance (or lack of it) as physical laws. The First Law, for example, is mean- 
ingless without the concept of “force,” a word Newton used in all three laws. In 
fact, standing alone, the First Law conveys a precise meaning only for zero force-, 
that is, a body remaining at rest or in uniform (i.e., unaccelerated, rectilinear) 
motion is subject to no force whatsoever. A body moving in this manner is 
termed a free body (or free particle) . The question of the frame of reference 
with respect to which the “uniform motion” is to be measured is discussed in the 
following section. 

In pointing out the lack of content in Newton’s First Law, Sir Arthur 
Eddington § observed, somewhat facetiously, that all the law actually says is that 
“every particle continues in its state of rest or uniform motion in a straight line 



*Truesdell (Tr68) points out that Leonhard Euler (1707—1783) clarified and developed the 
Newtonian concepts. Euler “put most of mechanics into its modern form” and “made mechanics 
simple and easy” (p. 106). 

fErnstMach (1838-1916) expressed his view in his famous book first published in 1883; E. Mach, Die 
Mechanic in ihrer Entwicklung historisch-kritisch dargestellt [The science of mechanics] (Prague, 1883). 
A translation of a later edition is available (Ma60). Interesting discussions are also given by 
R. B. Lindsay and H. Margeneau (Li36) and N. Feather (Fe59). 

^Enunciated in 1687 by Sir Isaac Newton (1642-1727) in his Philosophiae naturalis principia mathemal- 
ica [Mathematical principles of natural philosophy , normally called Principia \ (London, 1687). Previously, 
Galileo (1564-1642) generalized the results of his own mathematical experiments with statements 
equivalent to Newton’s First and Second Laws. But Galileo was unable to complete the description of 
dynamics because he did not appreciate the significance of what would become Newton’s Third 
Law — and therefore lacked a precise meaning of force. 

§Sir Arthur Eddington (Ed30, p. 124). 
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except insofar as it doesn’t.” This is hardly fair to Newton, who meant something 
very definite by his statement. But it does emphasize that the First Law by itself 
provides us with only a qualitative notion regarding “force.” 

The Second Law provides an explicit statement: Force is related to the time 
rate of change of momentum. Newton appropriately defined momentum (al- 
though he used the term quantity of motion) to be the product of mass and veloc- 
ity, such that 

p = mv (2.1) 



Therefore, Newton’s Second Law can be expressed as 



F = 



dp 

dt 




( 2 . 2 ) 



The definition of force becomes complete and precise only when “mass” is de- 
fined. Thus the First and Second Laws are not really “laws” in the usual sense; 
rather, they may be considered definitions. Because length, time, and mass are 
concepts normally already understood, we use Newton’s First and Second Laws 
as the operational definition of force. Newton’s Third Law, however, is indeed a 
laiv. It is a statement concerning the real physical world and contains all of the 
physics in Newton's laws of motion.* 

We must hasten to add, however, that the Third Law is not a general law of 
nature. The law does apply when the force exerted by one (point) object on an- 
other (point) object is directed along the line connecting the objects. Such 
forces are called central forces; the Third Law applies whether a central force is 
attractive or repulsive. Gravitational and electrostatic forces are central forces, 
so Newton’s laws can be used in problems involving these types of forces. 
Sometimes, elastic forces (which are actually macroscopic manifestations of mi- 
croscopic electrostatic forces) are central. For example, two point objects con- 
nected by a straight spring or elastic string are subject to forces that obey the 
Third Law. Any force that depends on the velocities of the interacting bodies is 
noncentral, and the Third Law may not apply. Velocity-dependent forces are 
characteristic of interactions that propagate with finite velocity. Thus the force 
between moving electric charges does not obey the Third Law, because the force 
propagates with the velocity of light. Even the gravitational force between mov- 
ing bodies is velocity dependent, but the effect is small and difficult to detect. 
The only observable effect is the precession of the perihelia of the inner planets 
(see Section 8.9). We will return to a discussion of Newton’s Third Law in 
Chapter 9. 

To demonstrate the significance of Newton’s Third Law, let us paraphrase it 
in the following way, which incorporates the appropriate definition of mass: 



*The reasoning presented here, viz., that the First and Second Laws are actually definitions and that 
the Third Law contains the physics, is not the only possible interpretation. Lindsay and Margenau 
(Li36), for example, present the first two Laws as physical laws and then derive the Third Law as a 
consequence. 
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III'. If two bodies constitute an ideal, isolated system, then the accelerations of these bodies 
are always in opposite directions, and the ratio of the magnitudes of the accelerations 
is constant. This constant ratio is the inverse ratio of the masses of the bodies. 



With this statement, we can give a practical definition of mass and therefore give 
precise meaning to the equations summarizing Newtonian dynamics. For two 
isolated bodies, 1 and 2, the Third Law states that 



F, = -F 2 

Using the definition of force as given by the Second Law, we have 

dp i rfp 2 

dt dt 



(2.3) 



(2.4a) 



or, with constant masses, 




and, because acceleration is the time derivative of velocity, 

wi(a,) = w 2 (- a 2 ) 



Hence, 



(2.4b) 



(2.4c) 



m 2 a i 

W] fl 2 



(2.5) 



where the negative sign indicates only that the two acceleration vectors are op- 
positely directed. Mass is taken to be a positive quantity. 

We can always select, say, m , as the unit mass. Then, by comparing the ratio 
of accelerations when m } is allowed to interact with any other body, we can de- 
termine the mass of the other body. To measure the accelerations, we must have 
appropriate clocks and measuring rods; also, we must choose a suitable coordi- 
nate system or reference frame. The question of a “suitable reference frame” is 
discussed in the next section. 

One of the more common methods of determining the mass of an object is 
by weighing — for example, by comparing its weight to that of a standard by 
means of a beam balance. This procedure makes use of the fact that in a gravita- 
tional field the weight of a body is just the gravitational force acting on the body; 
that is, Newton’s equation F = ma becomes W = mg, where g is the acceleration 
due to gravity. The validity of using this procedure rests on a fundamental as- 
sumption: that the mass m appearing in Newton’s equation and defined accord- 
ing to Statement III' is equal to the mass m that appears in the gravitational force 
equation. These two masses are called the inertial mass and gravitational mass, 
respectively. The definitions may be stated as follows: 

Inertial Mass: That mass determining the acceleration of a body under the action of a 
given force. 

Gravitational Mass: That mass determining the gravitational forces between a body 
and other bodies. 




52 



2 / NEWTONIAN MECHANICS— SINGLE PARTICLE 



Galileo was the First to test the equivalence of inertial and gravitational mass 
in his (perhaps apocryphal) experiment with falling weights at the Tower of Pisa. 

Newton also considered the problem and measured the periods of pendula of 
equal lengths but with bobs of different materials. Neither Newton nor Galileo \ 

found any difference, but the methods were quite crude.* In 1890 Eotvos* de- 
vised an ingenious method to test the equivalence of inertial and gravitational 
masses. Using two objects made of different materials, he compared the effect of 
the Earth’s gravitational force (i.e., the weight) with the effect of the inertial | 

force caused by the Earth’s rotation. The experiment involved a null method 
using a sensitive torsion balance and was therefore highly accurate. More recent 
experiments (notably those of Dicke*), using essentially the same method, have 
improved the accuracy, and we know now that inertial and gravitational mass are 
identical to within a few parts in 10 12 . This result is considerably important in the \ 

general theory of relativity. 8 The assertion of the exact equality of inertial and 1 

gravitational mass is termed the principle of equivalence. | 

Newton’s Third Law is stated in terms of two bodies that constitute an iso- 
lated system. It is impossible to achieve such an ideal condition; every body in the 
universe interacts with every other body, although the force of interaction may be 
far too weak to be of any practical importance if great distances are involved. 

Newton avoided the question of how to disentangle the desired effects from all 
the extraneous effects. But this practical difficulty only emphasizes the enormity 
of Newton’s assertion made in the Third Law. It is a tribute to the depth of his 
perception and physical insight that the conclusion, based on limited observa- 
tions, has successfully borne the test of experiment for 300 years. Only within the 
20th century did measurements of sufficient detail reveal certain discrepancies 
with the predictions of Newtonian theory. The pursuit of these details led to the 
development of relativity theory and quantum mechanics. 11 

Another interpretation of Newton’s Third Law is based on the concept of 
momentum. Rearranging Equation 2.4a gives 

^<Pi + P2) = o 
or 

Pi 4 P 2 — constant (2.6) \ 

The statement that momentum is conserved in the isolated interaction of two 
particles is a special case of the more general conservation of linear momen- 
tum. Physicists cherish general conservation laws, and the conservation of lin- 
ear momentum is believed always to be obeyed. Later we shall modify our defi- > 



*In Newton's experiment, he could have detected a difference of only one part in 10 s . 

fRoland von Eotvos (1848-1919), a Hungarian baron; his research in gravitational problems led to 

the development of a gravimeter, which was used in geological studies. 

JP. G. Roll, R. Krotkov, and R. H. Dicke, Ann. Phys. (N.Y.) 26, 442 (1964). See also Braginsky and 
Pavov, Sov. Phys.-JETP 34, 463 (1972). 

§See, for example, the discussions by P. G. Bergmann (Be46) and J. Weber (We61). Weber’s book 
also provides an analysis of the Eotvos experiment. 

IlSee also Section 2.8. 
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nition of momentum from Equation 2.1 for high velocities approaching the 
speed of light. 



2.3 Frames of Reference 

Newton realized that, for the laws of motion to have meaning, the motion of 
bodies must be measured relative to some reference frame. A reference frame is 
called an inertial frame if Newton’s laws are indeed valid in that frame; that is, if 
a body subject to no external force moves in a straight line with constant velocity 
(or remains at rest) , then the coordinate system establishing this fact is an iner- 
tial reference frame. This is a clear-cut operational definition and one that also 
follows from the general theory of relativity. 

If Newton’s laws are valid in one reference frame, then they are also valid in 
any reference frame in uniform motion (i.e., not accelerated) with respect to 
the first system.* This is a result of the fact that the equation F — mr involves the 
second time derivative of r: A change of coordinates involving a constant velocity 
does not influence the equation. This result is called Galilean invariance or the 
principle of Newtonian relativity. 

Relativity theory has shown us that the concepts of absolute rest and an ab- 
solute inertial reference frame are meaningless. Therefore, even though we con- 
ventionally adopt a reference frame described with respect to the “fixed” stars — 
and, indeed, in such a frame the Newtonian equations are valid to a high degree 
of accuracy — such a frame is, in fact, not an absolute inertial frame. We may, 
however, consider the “fixed” stars to define a reference frame that approxi- 
mates an “absolute” inertial frame to an extent quite sufficient for our present 
purposes. 

Although the fixed-star reference frame is a conveniently definable system 
and one suitable for many purposes, we must emphasize that the fundamental 
definition of an inertial frame makes no mention of stars, fixed or otherwise. If a 
body subject to no force moves with constant velocity in a certain coordinate sys- 
tem, that system is, by definition, an inertial frame. Because precisely describing 
the motion of a real physical object in the real physical world is normally diffi- 
cult, we usually resort to idealizations and approximations of varying degree; 
that is, we ordinarily neglect the lesser forces on a body if these forces do not sig- 
nificantly affect the body’s motion. 

If we wish to describe the motion of, say, a free particle and if we choose for 
this purpose some coordinate system in an inertial frame, then we require that 
the (vector) equation of motion of the particle be independent of the position of 
the origin of the coordinate system and independent of its orientation in space. 
We further require that time be homogeneous; that is, a free particle moving 
with a certain constant velocity in the coordinate system during a certain time 



*In Chapter 10, we discuss the modification of Newton’s equations that must be made if it is desired 
to describe the motion of a body with respect to a noninertial frame of reference, that is, a frame that 
is accelerated with respect to an inertial frame. 
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FIGURE 2-1 We choose to describe the path of a free particle moving along the path 
AC in a rectangular coordinate system whose origin moves in a circle. 
Such a system is not an inertial reference frame. 



interval must not, during a later time interval, be found to move with a different 
velocity. 

We can illustrate the importance of these properties by the following exam- 
ple. Consider, as in Figure 2-1, a free particle moving along a certain path AC. To 
describe the particle’s motion, let us choose a rectangular coordinate system 
whose origin moves in a circle, as shown. For simplicity, we let the orientation of 
the axes be fixed in space. The particle moves with a velocity relative to an in- 
ertial reference frame. If the coordinate system moves with a linear velocity v c 
when at the point B, and if \ c = \ p , then to an observer in the moving coordinate 
system the particle (at A) will appear to be at rest At some later time, however, 
when the particle is at C and the coordinate system is at D , the particle will ap- 
pear to accelerate with respect to the observer. We must, therefore, conclude 
that the rotating coordinate system does not quality' as an inertial reference 
frame. 

These observations are not sufficient to decide whether time is homoge- 
neous. To reach such a conclusion, repeated measurements must be made in 
identical situations at various times; identical results would indicate the homo- 
geneity of time. 

Newton’s equations do not describe the motion of bodies in noninertial sys- 
tems. We can derise a method to describe the motion of a particle by a rotating 
coordinate system, but, as we shall see in Chapter 10, the resulting equations con- 
tain several terms that do not appear in the simple Newtonian equation F = ma. 
For the moment, then, we restrict our attention to inertial reference frames to 
describe the dynamics of particles. 
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2.4 The Equation of Motion for a Particle 



Newton’s equation F = dp/dt can be expressed alternatively as 



F = 



— (mv) = m 
dt 



dv 

dt 



— mr 



( 2 . 7 ) 



if we assume that the mass m does not vary with time. This is a second-order dif- 
ferential equation that may be integrated to find r = r(!) if the function F is 
known. Specifying the initial values of r and r = v then allows us to evaluate the 
two arbitrary constants of integration. We then determine the motion of a parti- 
cle by the force function F and the initial values of position r and velocity v. 

The force F may be a function of any combination of position, velocity, and 
time and is generally denoted as F(r, v, t). For a given dynamic system, we nor- 
mally want to know r and v as a function of time. Solving Equation 2.7 will help 
us do this by solving for r. Applying Equation 2.7 to physical situations is an im- 
portant part of mechanics. 

In this chapter, we examine several examples in which the force function is 
known. We begin by looking at simple force functions (either constant or de- 
pendent on only one of r, v, and l) in only one spatial dimension as a refresher 
of earlier physics courses. It is important to form good habits in problem solving. 
Here are some useful problem-solving techniques. 

1. Make a sketch of the problem, indicating forces, velocities, and so forth. 

2. Write down the given quantities. 

3. Write down useful equations and what is to be determined. 

4. Strategy and the principles of physics must be used to manipulate the equa- 
tions to find the quantity sought. Algebraic manipulations as well as differ- 
entiation or integration is usually required. Sometimes numerical calcula- 
tions using a computer are the easiest, if not the only, method of solution. 

5. Finally, put in the actual values for the assumed variable names to determine 
the quantity sought. 

Let us first consider the problem of a block sliding on an inclined plane. Let 
the angle of the inclined plane be 6 and the mass of the block be 100 g. The 
sketch of the problem is shown in Figure 2-2a. 





FIGURE 2-2 Examples 2.1 and 2.2. 
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EXAMPLE, 2.1 



If a block slides without friction down a fixed, inclined plane with 0 — 30°, w'hat 
is the block’s acceleration? 

Solution. Two forces act on the block (see Figure 2-2a) : the gravitational force F 
and the plane’s normal force N pushing upward on the block (no friction in this 
example) . The block is constrained to be on the plane, and the only direction the 
block can move is the ^-direction, up and down the plane. We take the + x-direc- 
tion to be down the plane. The total force F net is constant; Equation 2.7 becomes 

F = F + N 

r net r g T 

and because F liet is the net resultant force acting on the block, 

Fnet = 
or 

r g + N = mr (2.8) 

This vector must be applied in two directions: x and y (perpendicular to x) . 
The component of force in the y-direction is zero, because no acceleration oc- 
curs in this direction. The force F is divided vectorially into its x- and y-compo- 
nents (dashed lines in Figure 2-2a). Equation 2.8 becomes 

y-direction 



-F g cos e + N = 0 



(2.9) 



x-direction 



F g sin 6 = mx (2.10) 

with the required result 

F e mp sin Q 

x = — sin 0 — — 2 sin 6 

m m 

x — g sin (30°) = - = 4.9 m/s 2 (2.11) 

Therefore the acceleration of the block is a constant. 

We can find the velocity of the block after it moves from rest a distance x 0 
down the plane by multiplying Equation 2.1 1 by 2x and integrating 



2xx = 2xg sin 6 



d, 

dt 



(x 2 ) = 2 g sin 6 



dx 

dt 



d(x 2 ) = 2g sin 6 



*0 



dx 



0 



'0 
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At t = 0, both x = x = 0, and, at t = < final , x = x 0 , and the velocity x = v 0 . 

Vo = 2 g sin 9 x 0 
v 0 = V2g sin 9 x 0 



EXAMPLE 2.2 



If the coefficient of static friction between the block and plane in the previous 
example is fx s = 0.4, at what angle 9 will the block start sliding if it is initially at 
rest? 

Solution. We need a new sketch to indicate the additional frictional force / (see 
Figure 2-2b). The static frictional force has the approximate maximum value 

fm ax - HsN ( 2 . 12 ) 

and Equation 2.7 becomes, in component form, 



y-direction 


—F g cos 9 + N — 0 


(2.13) 


x-direction 


— f + Fjj. sin 9 = mx 


(2.14) 



The static frictional force / will be some value / ^ / nax required to keep 5c = 0 
— that is, to keep the block at rest. However, as the angle 9 of the plane in- 
creases, eventually the static frictional force will be unable to keep the block at 
rest. At that angle 9', f becomes 

f{9 = O') = / niax = ix S N - fi s F g c os 9 

and 

mx = F g sin 9 - f max 

mx = F g sin 9 — fx s F g c os 9 (2.15) 

x = gfsin 9 — jx s cos 9) 

Just before the block starts to slide, the acceleration x = 0, so 

sin 9 — /x x cos 9 = 0 
tan 6 = /x s = 0.4 




After the block in the previous example begins to slide, the coefficient of ki- 
netic (sliding) friction becomes fx k = 0.3. Find the acceleration for the angle 
9 = 30°. 
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Solution. Similarly to Example 2.2, the kinetic friction becomes (approxi- 



mately) 




fk~ Vk N = Hk F g cos 8 


(2.16) 


and 




mx = Fg sin 6 — f k = mg (sin 9 — fx k cos 9) 


(2.17) 




x = g (sin 9 — n h cos 0) = 0.24 g 


(2.18) 



Generally, the force of static friction (f max = /jl s N ) is greater than that of 
kinetic friction (J k = fx k N). This can be observed in a simple experiment. If we 
lower the angle 9 below 16.7°, we find that x < 0, and the block eventually 
stops. If we raise the block back up above 6 = 16.7°, we find that the block does 
not start sliding again until 6 > 22° (Example 2.2). The static friction deter- 
mines when it starts moving again. There is not a discontinuous acceleration as 
the block starts moving, because of the difference between /r 5 and jx k . For small 
speeds, the coefficient of friction changes rather quickly from yu ( to fx k . 

The subject of friction is still an interesting and important area of research. 
There are still surprises. For example, even though we calculate the absolute 
value of the frictional force as /= fiN, research has shown that the frictional 
force is directly proportional, not to the load, but to the microscopic area of 
contact between the two objects (as opposed to the apparent contact area). We 
use fiN as an approximation because, as N increases, so does the actual contact 
area on a microscopic level. For hundreds of years before the 1940s, it was ac- 
cepted that the load — and not the area — was directly responsible. We also be- 
lieve that the static frictional force is larger than that of kinetic friction because 
the bonding of atoms between the two objects does not have as much time to 
develop in kinetic motion. 



Effects of Retarding Forces 

We should emphasize that the force F in Equation 2.7 is not necessarily constant, 
and indeed, it may consist of several distinct parts, as seen in the previous exam- 
ples. For example, if a particle falls in a constant gravitational field, the gravita- 
tional force is F g = mg, where g is the acceleration of gravity. If, in addition, a 
retarding force F, exists that is some function of the instantaneous speed, then 
the total force is 



F = f + p 

M. R g I A r 

~= mg + F,(u) 



(2.19) 



It is frequently sufficient to consider that F r (u) is simply proportional to some 
power of the speed. In general, real retarding forces are more complicated, but 
the power-law approximation is useful in many instances in which the speed 
does not vary greatly. Even more to the point, if F r oc v n , then the equation of 
motion can usually be integrated directly, whereas, if the true velocity depend- 
ence were used, numerical integration would probably be necessary. With the 
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power-law approximation, we can then write 

F — mg — mkv n - (2.20) 

where k is a positive constant that specifies the strength of the retarding force 
and where \/v is a unit vector in the direction of v. Experimentally, we find 
that, for a relatively small object moving in air, n = 1 for velocities less than 
about 24 m/s (~ 80 ft/s). For higher velocities but below the velocity of sound 
(~ 330 m/s or 1,100 ft/s), the retarding force is approximately proportional to 
the square of the velocity.* For simplicity, the v 2 dependence is usually taken 
for speeds up to the speed of sound. 

The effect of air resistance is important for a ping-pong ball smashed to an 
opponent, a high-flying softball hit deep to the outfield, a golfer’s chip shot, and 
a mortar shell lofted against an enemy. Extensive tabulations have been made 
for military ballistics of projectiles of various sorts for the velocity as a function of 
flight time. There are several forces on an actual projectile in flight. The air re- 
sistance force is called the drag W and is opposite to the projectile’s velocity as 
shown in Figure 2-3a. The velocity v is normally not along the symmetry axis of 
the shell. The component of force acting perpendicular to the drag is called the 
lift L n . There may also be various other forces due to the projectile’s spin and os- 
cillation, and a calculation of a projectile’s ballistic trajectory is quite complex. 
The Prandtl expression for the air resistance^ is 

1 2 

W=~c w pAv* (2.21) 

where c w is the dimensionless drag coefficient, p is the air density, v is the veloc- 
ity, and A is the cross-sectional area of the object (projectile) measured perpen- 
dicularly to the velocity. In Figure 2-3b, we plot some typical values for c w , and in 
Figures 2-3c and d we display the calculated air resistance Wusing Equation 2.21 
for a projectile diameter of 10 cm and using the values of c w shown. The air re- 
sistance increases dramatically near the speed of sound (Mach number M = 
speed/speed of sound). Below speeds of about 400 m/s it is evident that an 
equation of at least second degree is necessary to describe the resistive force. For 
higher speeds, the retarding force varies approximately linearly with speed. 

Several examples of the motion of a particle subjected to various forces are 
given below. These examples are particularly good to begin computer calcula- 
tions using any of the available commercial math programs and spreadsheets or 
for the students to write their own programs. The computer results, especially 
the plots, can often be compared with the analytical results presented here. 
Some of the figures shown in this section were produced using a computer, and 



*The motion of a particle in a medium in which there is a resisting force proportional to the speed 
or to the square of the speed (or to a linear combination of the two) was examined by Newton in his 
Principia (1687). The extension to any power of the speed was made by Johann Bernoulli in 1711. 
The term Stokes' laxo of resistance is sometimes applied to a resisting force proportional to the speed; 
Newton’s law of resistance is a retarding force proportional to the square of the speed. 
tSee the article by E. Melchior and M. Reuschel in Handbook on Weaponry (Rh82, p. 137). 
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Velocity (m/s) Velocity (m/s) 



(c) 



<d) 



FIGURE 2-3 (a) Aerodynamic forces acting on projectile. W is the drag (air resistive 

force) and is opposite the velocity of the projectile v. Notice that v may 
be at an angle a from the symmetry axis of projectile. The component 
of force acting perpendicular to the drag is called the lift L a . The point 
D is the center of pressure. Finally, the gravitational force F g acts down. 
If the center of pressure is not at the projectile's center of mass, there is 
also a torque about the center of mass, (b) The drag coefficient c w , 
from the Rheinmetall resistance law (Rh82), is plotted versus the Mach 
number M. Notice the large change near the speed of sound where 
M = 1. (c) The air resistive force W (drag) is shown as a function of 
velocity for a projectile diameter of 10 cm. Notice the inflection near 
the speed of sound, (d) Same as (c) for higher velocities. 



several end-of-chapter problems are meant to develop the student’s computer 
experience if so desired by the instructor or student. 



EXAMPLE 2.4 



As the simplest example of the resisted motion of a particle, find the displace- 
ment and velocity of horizontal motion in a medium in which the retarding 
force is proportional to the velocity. 

Solution. A sketch of the problem is shown in Figure 2-4. The Newtonian equa- 
tion F — ma provides us with the equation of motion: 
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Resisting force F = kmv 



FIGURE 2-4 Example 2.4. 



x-direction 



dv 

ma — m — = — kmv 
dt 



( 2 . 22 ) 



where kmv is the magnitude of the resisting force ( k = constant). We are not 
implying by this form that the retarding force depends on the mass m; this form 
simply makes the math easier. Then 




(2.23) 



Inn = — kt + Q 

The integration constant in Equation 2.23 can be evaluated if we prescribe the 
initial condition v{t = 0) = v 0 . The Q = In v 0 , and 

v = v 0 e~ kt ' (2.24) 



We can integrate this equation to obtain the displacement x as a function of 
time: 



v 



dx 



dt 



= v 0 e 



x 



I'o 



e~ kt dt = -—e~ kt + C 2 

« 



(2.25a) 



The initial condition x(t = 0) = 0 implies C 2 = v Q /k. Therefore 

* = “ e ~ kt ) (2.25b) 

fe 

This result shows that x asymptotically approaches the value v 0 /k as t—>oo. 

We can also obtain the velocity as a function of displacement by writing 

dv _ dv dt _ dv 1 
dx dt dx dt v 



so that 



dv _ dv 
dx dt 



— kv 



dv 

dx 



-k 



or 
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from which we find, by using the same initial conditions, 

v — V[) — kx 

Therefore, the velocity decreases linearly with displacement. 



(2.26) 




Find the displacement and velocity of a particle undergoing vertical motion in a 
medium having a retarding force proportional to the velocity. 

Solution. Let us consider that the particle is falling downward with an initial 
velocity Vq from a height hin a constant gravitational field (Figure 2-5). The 
equation of motion is 



z-direction 



F = m — = — rag — kmv 
dt 



(2.27) 



where — kmv represents a positive upward force since we take z and v = z to be 
positive upward, and the motion is downward — that is, v < 0 , so that —kmv > 0. 
From Equation 2.27, we have 

dv 

(2.28) 

kv + g 

Integrating Equation 2.28 and setting v(t = 0) = v 0 , we have (noting that v 0 < 0) 



\n(kv + g) = —t + c 
kv + g = e kt+kc 



v = — = -- + 
dt k 



g , *«/ 0 + g 



(2.29) 




FIGURE 2-5 Example 2.5. 
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FIGURE 2-6 Results for Example 2.5 indicating the downward speeds for various 
initial speeds v 0 as they approach the terminal velocity. 



Integrating once more and evaluating the constant by setting z(t — 0) = h, we 
find 



fit kVa + fi 

h-\+ ■ g (1 — e~ kt ) ■ (2.30) 

k k z 

Equation 2.29 shows that as the time becomes very long, the velocity ap- 
proaches the limiting value —g/k\ this is called the terminal velocity, v t . 
Equation 2.27 yields the same result, because the force will vanish — and hence 
no further acceleration will occur — when v = ~g/k. If the initial velocity ex- 
ceeds the terminal velocity in magnitude, then the body immediately begins to 
slow down and v approaches the terminal speed from the opposite direction. 
Figure 2-6 illustrates these results for the downward speeds (positive values). 



EXAMPLE 2.6 



Next, we treat projectile motion in two dimensions, first without considering air 
resistance. Let the muzzle velocity of the projectile be u 0 and the angle of eleva- 
tion be 6 (Figure 2-7). Calculate the projectile’s displacement, velocity, and range. 

Solution. Using F = mg, the force components become 
x-dirextion 



0 = nix 



(2.31a) 



y-direction 



— mg = my 



(2.31b) 
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FIGURE 2-7 Example 2.6. 



Neglect the height of the gun, and assume x = y = 0 at t = 0. Then 

x — 0 

x = v 0 cos 8 

x = v 0 t cos 8 (2.32) 

and 



y = -g 

y = ~gt + v 0 sin 8 

— gt 2 

y = — 1- v 0 t sin 8 (2.33) 

The speed and total displacement as functions of time are found to be 

v = \/ x 2 + y 2 = {v\ + g 2 t 2 — gt sin d ) l/ 2 (2.34) 

and 

/ / g-t 2 V/ 2 

r — v x 2 + y 2 = ( v 2 t 2 -I — r> o gt 3 sin0j (2.35) 



We can find the range by determining the value of xwhen the projectile falls 
back to ground, that is, when y = 0. 



y = 1 



+ v 0 sin 0 I = 0 



One value of y — 0 occurs for t — 0 and the other one for t — T. 



+ v 0 sin 0—0 



2v 0 sin 0 
g 



(2.37) 
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The range R is found from 



x(t = T) ~ range 




sin 9 cos 6 



R = range 



— sin 29 
g 



(2.38) 

(2.39) 



Notice that the maximum range occurs for 6 = 45°. 

Let us use some actual numbers in these calculations. The Germans used a 
long-range gun named Big Bertha in World War I to bombard Paris. Its muzzle 
velocity was 1,450 m/s. Find its predicted range, maximum projectile height, 
and projectile time of flight if 6 — 55°. We have u 0 = 1450 m/s and 6 = 55°, so 
the range (from Equation 2.39) becomes 



R = 



(1450 m/s) 2 
9.8 m/s 2 



[sin(110 0 )] - 202 km 



Big Bertha’s actual range was 120 km. The difference is a result of the real 
effect of air resistance. 

To find the maximum predicted height, we need to calculated y for the 
time 7/2 where Tis the projectile time of flight: 



= (2) (1450 m/s) (sin 55°) 

9.8 m/s 2 S 

/ T\ ~ gT 2 v 0 T . 

= ^ I = — g- + — Sin 6 

— (9.8 m/s) (242 s) 2 i (1450 m/s) (242 s) sin(55°) 
4 - 



8 



= 72 km 



EXAMPLE 2.7 



Next, we add the effect of air resistance to the motion of the projectile in the 
previous example. Calculate the decrease in range under the assumption that 
the force caused by air resistance is directly proportional to the projectile’s 
velocity. 

Solution. The initial conditions are the same as in the previous example. 

x(t = 0) = 0 = y{t = 0) 1 

x(t = 0) — v 0 cos 6 = U > (2.40) 

y{t— 0) = u 0 sin 9 ss V J 

However, the equations of motion, Equation 2.31, become 

mx — —kmx 
my = —kmy ~ mg 



(2.41) 

(2.42) 
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Equation 2.41 is exactly that used in Example 2.4. The solution is therefore 



x — 



U 

k 



(1 - e~ kt ) 



( 2 . 43 ) 



Similarly, Equation 2.42 is the same as the equation of the motion in Example 
2.5. We can use the solution found in that example by letting h — 0. (The fact 
that we considered the particle to be projected downward in Example 2.5 is of no 
consequence. The sign of the initial velocity automatically takes this into ac- 
count.) Therefore 



>’ 




kV+ g 
k 2 



(1 - e ~ kt ) 



( 2 . 44 ) 



The trajectory is shown in Figure 2-8 for several values of the retarding force 
constant k for a given projectile flight. 

The range R', which is the range including air resistance, can be found as 
previously by calculating the time T required for the entire trajectory and then 
substituting this value into Equation 2.43 for x. The time T is found as previ- 
ously by finding i — T when y = 0. From Equation 2.44, we find 

kV+g 

T= - e~ kT ) ( 2 . 45 ) 

gk 

This is a transcendental equation, and therefore we cannot obtain an analytic 
expression for T. Nonetheless, we still have powerful methods to use to solve 



y 




FIGURE 2-8 The calculated trajectories of a particle in air resistance (F res = —kmv) 

for various values of k (in units of s -1 ). The calculations were performed 
for values of 6 = 60° and v 0 — 600 m/s. The values of y (Equation 2.44) 
are plotted versus x (Equation 2.43). 
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such problems. We present two of them here: (1 ) a perturbation method to find an 
approximate solution, and (2) a numerical method, which can normally be as ac- 
curate as desired. We will compare the results. 



Perturbation Method To use the perturbation method, we find an expansion pa- 
rameter or coupling constant that is normally small. In the present case, this param- 
eter is the retarding force constant k, because we have already solved the present 
problem with k = 0, and now we would like to turn on the retarding force, but 
let k be small. We therefore expand the exponential term of Equation 2.45 (see 
Equation D.34 of Appendix D) in a power series with the intention of keeping 
only the lowest terms of k n , where h is our expansion parameter. 



T = 



kV+ g 
g k 



1 1 

kT - - k 2 T 2 + -k‘ 3 T $ - 

, 2 6 



(2.46) 



If we keep only terms in the expansion through A 3 , this equation can be re- 
arranged to yield 



T = 



2 V/g 1 9 

— 1 — kT 2 

1 + kV/g 3 



(2.47) 



We now have the expansion parameter k in the denominator of the first term on 
the right-hand side of this equation. We need to expand this term in a power se- 
ries (Taylor series, see Equation D.8 of Appendix D): 



YVWg = 1 “ {kv/s) + {kV/g) 2 ~ (2 ‘ 48) 

where we have kept only terms through k 2 , because we only have terms through 
k in Equation 2.47. If we insert this expansion of Equation 2.48 into the first 
term on the right-hand side of Equation 2.47 and keep only the terms in k to first 
order, we have 



T = 






k + 0{ k 2 ) 



(2.49) 



where we choose to neglect 0{k 2 ), the terms of order kr and higher. In the limit 
k — > 0 (no air resistance), Equation 2.49 gives us the same result as in the previ- 
ous example: 



T(k = 0) 



To 



2V 2v 0 sin 6 
g g 



Therefore, if k is small (but nonvanishing), the flight time will be approximately 
equal to T 0 . If we then use this approximate value for T — T 0 in the right-hand 
side of Equation 2.49, we have 



T = 





(2.50) 



which is the desired approximate expression for the flight time. 
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Next, we write the equation for x (Equation 2.43) in expanded form: 

x = — — — k 2 / 2 + - k*t $ — • • ■ ^ ( 2 . 51 ) 

Because x{t = T) = R', we have approximately for the range 

R' s u(t- ~kT^j ( 2 . 52 ) 



where again we keep terms only through the first order of k. We can now evalu- 
ate this expression by using the value of T from Equation 2.50. If we retain only 
terms linear in k, we find 



g V $g) 



( 2 . 53 ) 



The quantity 26 r V/gcan now be written (using Equations 2.40) as 

2 UV 2 vl vl 

= sin 6 cos 6 = — sin 2 8 = R (2.54) 

g g g 



which will be recognized as the range R of the projectile when air resistance is 
neglected. Therefore 




( 2 . 55 ) 



Over what range of values for k would we expect our perturbation method to be 
correct? If we look at the expansion in Equation 2.48, we see that the expansion 
will not converge unless kV/g < 1 or k < g/V, and in fact, we would like k <5C 
g/V= g/(v o sin 6). 



Numerical Method Equation 2.45 can be solved numerically using a computer 
by a variety of methods. We set up a loop to solve the equation for T for many 
values of k up to 0.08 s _1 : T t (k t ). These values of T t and k t are inserted into 
Equation 2.43 to find the range R' t , which is displayed in Figure 2-9. The range 
drops rapidly for increased air resistance, just as one would expect, but it does 
not display the linear dependence suggested by the perturbation method solu- 
tion of Equation 2.55. 

For the projectile motion described in Figures 2-8 and 2-9, the linear ap- 
proximation is inaccurate for k values as low as 0.01 s” 1 and incorrectly shows 
the range is zero for all values of k larger than 0.014 s _1 . This disagreement with 
the perturbation method is not surprising because the linear result for the range 
R' was dependent on k <SC g/(v 0 sin 6) = 0.02 s ~ 1 , which is hardly true for even 
k = 0.01 s -1 . The agreement should be adequate for k = 0.005 s _l . The results 
shown in Figure 2-8 indicate that for values of k > 0.005 s -1 , the drag can hardly 
be considered a perturbation. In fact, for k > 0.01 s -1 the drag becomes the 
dominant factor in the projectile motion. 
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Retarding force constant, k (s 1 ) 

FIGURE 2-9 The range values calculated approximately and numerically for the 
projectile data given in Figure 2-8 are plotted as a function of the 
retarding force constant k. 



The previous example indicates how complicated the real world can be. In that 
example, we still had to make assumptions that were nonphysical — in assuming, 
for example, that the retarding force is ahvays linearly proportional to the veloc- 
ity. Even our numerical calculation is not accurate, because Figure 2-3 shows us 
that a better assumption would be to include a xr retarding term as well. 

Adding such a term would not be difficult with the numerical calculation, and 
we shall do a similar calculation in the next example. We have included the au- 
thor’s Mathcad file that produced Figures 2-8 and 2-9 in Appendix H for those 
students who might want to reproduce the calculation. We emphasize that there 
are many ways to perform numerical calculations with computers, and the stu- 
dent will probably want to become proficient with several. 



EXAMPLE 2.8 



Use the data shown in Figure 2-3 to calculate the trajectory for an actual pro- 
jectile. Assume a muzzle velocity of 600 m/s, gun elevation of 45°, and a pro- 
jectile mass of 30 kg. Plot the height y versus the horizontal distance x and 
plot y, x , and y versus time both with and without air resistance. Include only 
the air resistance and gravity, and ignore other possible forces such as the 
lift. 

Solution . First, we make a table of retarding force versus velocity by reading 
Figure 2-3. Read the force every 50 m/s for Figure 2-3c and every 100 m/s for 
Figure 2-3d. We can then use a straight line interpolation between the tabular 
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values. We use the coordinate system shown in Figure 2-7. The equations of 
motion become 



X 

m 


(2.56) 


F v 




(2.57) 



where F x and F y are the retarding forces. Assume gis constant. F x will always be a 
positive number, but F y > 0 for the projectile going up, and F y < 0 for the pro- 
jectile coming back down. Let 6 be the projectile’s elevation angle from the hor- 
izontal at any instant. 


v — VV 2 + y 1 


(2.58) 


y 

tan 9 = - 

X 


(2.59) 


F x = F cos 6 


(2.60) 


F y = Fsin 6 


(2.61) 


We can calculate F x and F y at any instant by knowing xand y. Over a small time 
interval, the next x and y can be calculated. 


. f ‘ , 

x — x at + v ( ) cos 6 

JO 


(2.62) 


P 

y — y dt + v 0 sin 6 
Jo 


(2.63) 


X 

II 

Sr 


(2.64) 


V; 

II 

Vi. 

sx 


(2.65) 



'o 



We wrote a short computer program to contain our table for the retarding 
forces and to perform the calculations for x, y, x, and y as a function of time. We 
must perform the integrals by summations over small time intervals, because 
the forces are time dependent. Figure 2-10 shows the results. 

Notice the large difference that the air resistance makes. In Figure 2- 10a, 
the horizontal distance (range) that the projectile travels is about 16 km com- 
pared to almost 37 km with no air resistance. Our calculation ignored the fact 
that the air density depends on the altitude. If we take account of the decrease 
in the air density with altitude, we obtain the third curve with a range of 18 km 
shown in Figure 2-1 0a. If we also included the lift, the range would be still 
greater. Notice that the change in velocities in Figures 2-10c and 2-10d mirror 
the air resistive force of Figure 2-3. The speeds decrease rapidly until the speed 
reaches the speed of sound, and then the rate of change of the speeds levels off 
somewhat. 
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(C) (d) 

FIGURE 2-10 The results of Example 2.8. The solid lines are the results if no air resis- 
tance is included, whereas the dashed lines include the results of adding 
the air resistive force. In (a) we also include the effect of the air density 
dependence, which becomes smaller as the projecule rises higher. 



This concludes our subsection on the effects of retarding forces. Much more 
could be done to include realistic effects, but the method is clear. Normally, one ef- 
fect is added at a time, and the results are analyzed before another effect is added. 

Other Examples of Dynamics 

We conclude this section with two additional standard examples of dynamical 
particle-like behavior. 



EXAMPLE 2.9 



Atwood’s machine consists of a smooth pulley with two masses suspended from 
a light string at each end (Figure 2-11). Find the acceleration of the masses and 
the tension of the string (a) when the pulley center is at rest and (b) when the 
pulley is descending in an elevator with constant acceleration a. 
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(a) (b) 

FIGURE 2-11 Example 2.9; Atwood’s machine. 



Solution. We neglect the mass of the string and assume that the pulley is 
smooth — that is, no friction on the string. The tension T must be the same 
throughout the string. The equations of: motion become, for each mass, for 
case (a), 

m x Xi — m } g — T (2.66) 

m 2 x = m 2 g ~ T (2.67) 



Notice again the advantage of the force concept: We need only identify the 
forces acting on each mass. The tension Tis the same in both equations. If 
the string is inextensible, then x 2 — ~ 3q, and Equations 2.66 and 2.67 may be 
combined 



mjX, - m ; g — ( m 2 g — m 2 x 2 ) 
= g- ( m 2 g+ m 2 x l ) 

Rearranging, 

g(m } - m 2 ) 

x , = = -* 2 

wt] + m 2 



( 2 . 68 ) 



If wi] > m 2 , then jq > 0, and x% < 0. The tension can be obtained from 
Equations 2.68 and 2.66: 



T = 
T = 



m x g — m x x x 



m \g ~ m ig 



{nil ~ m 2 ) 

m l + m 2 



2m 1 m 2 g 
raj + m 2 



(2.69) 
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For case (b), in which the pulley is in an elevator, the coordinate system 
with origins at the pulley center is no longer an inertial system. We need an in- 
ertial system with the origin at the top of the elevator shaft (Figure 2-1 lb). The 
equations of motion in the inertial system ( x" — x[ + x,, x” = x 2 + x 2 ) are 



so 



where x[ 

ing T: 



x 2 = 





•• // _ 
mjXj — 


m 


i(«i + so 




jf 

m 2 x 2 


TTh (x 2 + x 2 ) 


rn jXj 


= m lg ~ 


T 


»• t 

— m x X} — 


m 2 x 2 


= m 2 g~ 


T 


•• / 

— m 2 x < 2 — 


a. We have x 2 


— 


— x\, so we 



= ™\g ~ T 
= m 2 g~ 1 



m \ (g - a) - Tl 
m 2 (g~ aj — T) 



( 2 . 70 ) 



solve for xj as before by eliminat- 



“ x 2 - (g~ 



a) 



(m x — m 2 ) 
m [ + m 2 



( 2 . 71 ) 



and 



T = 



2m i m 2 (g - a) 
mi + m 2 



( 2 . 72 ) 



Notice that the results for the acceleration and tension are just as if the acceler- 
ation of gravity were reduced by the amount of the elevator acceleration a. 

The change for an ascending elevator should be obvious. 



EXAMPLE 2.10 



In our last example in this lengthy review of the equations of motion for a parti- 
cle, let us examine particle motion in an electromagnetic field. Consider a 
charged particle entering a region of uniform magnetic field B — for example, 
the earth’s field — as shown in Figure 2-12. Determine its subsequent motion. 

Solution. Choose a Cartesian coordinate system with its y-axis parallel to the 
magnetic field. If q is the charge on the particle, v its velocity, a its acceleration, 
and B the earth’s magnetic field, then 

v — xi + yj + zk 
a = xi + yj + zk 

B = B 0 j 

The magnetic force F = qv x B = ma, so 

m(x i + yj + zk) = ^(xi + yj + zk) X B 0 j — qB 0 (x k — zi) 
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z 




FIGURE 2~I2 Example 2.10; a moving particle enters a region of magnetic field. 



Equating like vector components gives 



mx ~ - qB^z j 

my = 0 > (2.73) 

m'i ~ qB 0 x ) 

Integrating the second of these equations, my = 0, yields 

5 ' = h 

where yo is a constant and is the initial value of y. Integrating a second time gives 

y = yot + yo 

where y 0 is also a constant. 

To integrate the first and last equations of Equation 2.73, let a = so 

that 



•« • 

x — ~az 

z = ax 



(2.74) 



These coupled, simultaneous differential equations can be easily uncoupled by 
differentiating one and substituting it into the other, giving 



z ~ ax ~ ~a 2 z 
x = — az ~ — a x 



so that 



z - — a‘z 



x 



—a 2 x 



(2.75) 
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Both of these differential equations have the same form of solution. Using the 
technique of Example C.2 of Appendix C, we have 

x — A cos at + B sin at + x 0 
z = A' cos at + B' sin at + z 0 

where A, A', B, B', x 0 , and z 0 are constants of integration that are determined by 
the particle’s initial position and velocity and by the equations of motion, 
Equation 2.74. These solutions can be rewritten 

(x — x 0 ) = A cos at + B sin at 1 

(y ~ y 0 ) = Jo 1 / ( 2 . 76 ) 

(z — z 0 ) = A' cos at + B' sin at j 

The x- and z-coordinates are connected by Equation 2.74, so substituting 
Equations 2.76 into the first equation of Equation 2.74 gives 

—a 2 A cos at — a 2 B sin at = —a(~aA' sin at + aB' cos at) ( 2 . 77 ) 

Because Equation 2.77 is valid for all t, in particular t = 0 and t ~ tt/2 a, 
Equation 2.77 yields 



so that 



— a 2 A = - a 2 B ' 
A = B' 



and 



gives 



— a 2 B — a 2 A' 



B= -A' 



We now have 



(x — x 0 ) = A cos at + B sin at 

(y ~ yo) = >’o t 

(z — z 0 ) = —B cos at A A sin at 



( 2 . 78 ) 



If at t — 0, z = z 0 and x = 0, then from Equation 2.78, differentiating and set- 
ting l = 0 gives 



aB = 0 



and 



aA — z 0 
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SO 



Finally, 



( x — x 0 ) = — cos at 
a 

(y - yo) = yot 

Z o 

(2 — z 0 ) = ~ sin at 



I 


( z (l mf 


| 


/ qB 0 t\ 


1 

X 

0 

1! 


\<j B oJ 


cos 


\m J 


O 

•**\ 

II 

0 

1 






II 

0 

1 


( z 0 m\ 


i • / 


'qBot\ 


\<l B 0 / 


Isinl 


v m ) 



( 2 . 79 ) 



These are the parametric equations of a circular helix of radius z 0 m/qB 0 . Thus, 
the faster the particle enters the field or the greater its mass, the larger the 
radius of the helix. And the greater the charge on the particle or the stronger 
the magnetic field, the tighter the helix. Notice also how the charged particle is 
captured by the magnetic field — just drifting along the field direction. In this 
example, the particle had no initial component of its velocity along the x-axis, 
but even if it had it would not drift along this axis (see Problem 2-31). Finally, 
notice that the magnetic force on the particle always acts perpendicular to its 
velocity and hence cannot speed it up. Equation 2.79 verifies this fact. 

The earth’s magnetic field is not as simple as the uniform field of this exam- 
ple. Nevertheless, this example gives some insight into one of the mechanisms by 
which the earth’s magnetic field traps low-energy cosmic rays and the solar wind to 
create the Van Allen belts. 



2.5 Conservation Theorems 

We now turn to a detailed discussion of the Newtonian mechanics of a single 
particle and derive the important theorems regarding conserved quantities. We 
must emphasize that we are not proving the conservation of the various quanti- 
ties. W T e are merely deriving the consequences of Newton’s laws of dynamics. 
These implicadons must be put to the test of experiment, and their verification 
then supplies a measure of confirmation of the original dynamical laws. The fact 
that these conservation theorems have indeed been found to be valid in many 
instances furnishes an important part of the proof for the correctness of 
Newton’s laws, at least in classical physics. 

The first of the conservation theorems concerns the linear momentum of a 
particle. If the particle is free, that is, if the particle encounters no force, then 
Equation 2.2 becomes simply p = 0. Therefore, p is a vector constant in time, 
and the first conservation theorem becomes 
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I. The total linear momentum p of a particle is conserved when the total force on it is 
zero. 

Note that this result is derived from a vector equation, p = 0, and therefore 
applies for each component of the linear momentum. To state the result in 
other terms, we let s be some constant vector such that F • s = 0, independent of 
time. Then 



p • s = F • s — 0 

or, integrating with respect to time, 

p * s = constant (2.80) 



which states that the componen t of linear momentum in a direction in which the force 
vanishes is constant in time. 

The angular momentum L of a particle with respect to an origin from which 
the position vector r is measured is defined to be 



L = r x p 



(2.81) 



The torque or moment of force N with respect to the same origin is defined 
to be 



N = r X F 



(2.82) 



where r is the position vector from the origin to the point where the force F is 
applied. Because F = mv for the particle, the torque becomes 

N = r x mv = rxp 



Now 



but 



d 

L = x P) = ( r x P> + ( r x P) 

r X p = r X mv = m(r X r) = 0 



so 



L = rXp = N 



(2.83) 



If no torques act on a particle (i.e., if N = 0), then L = 0 and L is a vector con- 
stant in time. The second important conservation theorem is 



II. The angular momentum of a particle subject to no torque is conserved. 

We remind the student that a judicious choice of the origin of a coordinate 
system will often allow a problem to be solved much more easily than a poor 
choice. For example, the torque will be zero in coordinate systems centered 
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along the resultant line of force. The angular momentum will be conserved in 
this case. 

If work is done on a particle by a force F in transforming the particle from 
Condition 1 to Condition 2, then this work is defined to be 

(2.84) 

If F is the net resultant force acting on the particle, 

dv dr dx 

F • dr — m — dt ~ m — • xdt 

dt dt dt 

= f |( vv )dt = f 2 /v*)d< = d(±mvj (2.85) 

The integrand in Equation 2.84 is thus an exact differential, and the work done 
by the total force F acting on a particle is equal to its change in kinetic energy: 

WJa = = ^m(vl - vf) = T, — 7\ (2.86) 

where T = ~mv 2 is the kinetic energy of the particle. If 7\ > then W 12 < 0, 
and the particle has done work with a resulting decrease in kinetic energy. It is 
important to realize that the force F leading to Equation 2.85 is the total (i.e., 
net resultant) force on the particle. 

Let us now examine the integral appearing in Equation 2.84 from a differ- 
ent standpoint. In many physical problems, the force F has the property that the 
work required to move a particle from one position to another without any 
change in kinetic energy depends only on the original and final positions and 
not on the exact path taken by the particle. For example, assume the work done 
to move the particle from point 1 in Figure 2-13 to point 2 is independent of the 
actual paths a, b, or c taken. This property is exhibited, for example, by a con- 
stant gravitational force field. Thus, if a particle of mass m is raised through a 
height h (by any path) , then an amount of work mgh has been done on the parti- 
cle, and the particle can do an equal amount of work in returning to its original 
position. This capacity to do work is called the potential energy of the particle. 

We may define the potential energy of a particle in terms of the work (done 
by the force F) required to transport the particle from a point 1 to a point 2 
(with no net change in kinetic energy): 

-2 

F-dr=Ui~U« (2.87) 

. i 

The work done in moving the particle is thus simply the difference in the poten- 
tial energy U at the two points. For example, if we lift a suitcase from position 1 
on the ground to position 2 in a car trunk, we as the external agent are doing 
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FIGURE 2-13 For some forces (identified later as conservative), the work done by the 
force to move a particle from one position 1 to another position 2 is 
independent of the path {a, b, or c). 



work against the force of gravity. Let the force F in Equation 2.87 be the gravita- 
tional force, and in raising the suitcase, F • dr becomes negative. The result of 
the integration in Equation 2.87 is that U x — U 2 is negative, so that the potential 
energy at position 2 in the car’s trunk is greater than that at position 1 on the 
ground. The change in potential energy U 2 — U x is the negative of the work 
done by the gravitational force, as can be seen by multiplying both sides of 
Equation 2.87 by —1. As the external agent, we do positive work (against gravity) 
to raise the potential energy' of the suitcase. 

Equation 2.87 can be reproduced* if we write F as the gradient of the scalar 
function U: 



Then 



F = -grad U = ~VU 



( 2 . 88 ) 



[ F* dr 



Ji 



n 



(V ( 7 ) * dr — - 



f 2 



dU — L\ - U 2 



( 2 . 89 ) 



In most systems of interest, the potential energy is a function of position 
and, possibly, time: U = U( r) or U — U( r, t). We do not consider cases in which 
the potential energy is a function of the velocity.* 

It is important, to realize that the potential energy is defined only to within an 
additive constant; that is, the force defined by — V U is no different from that de- 
fined by —V(U + constant). Potential energy therefore has no absolute meaning; 
only differences of potential energy are physically meaningful (as in Equation 2.87). 



*The necessary and sufficient condition that permits a vector function to be represented by the gra- 
dient of a scalar function is that the curl of the vector function vanishes identically. 
•(Velocity-dependent potentials are neccessary in certain situations, c.g., in electromagnetism (the 
so-called Lienard-Wiechert potentials). 
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If we choose a certain inertial frame of reference to describe a mechanical 
process, the laws of motion are the same as in any other reference frame in uni- 
form motion relative to the original frame. The velocity of a particle is in general 
different depending on which inertial reference frame we chose as the basis for 
describing the motion. We therefore find that it is impossible to ascribe an ab- 
solute kinetic energy to a particle in much the same way that it is impossible to 
assign any absolute meaning to potential energy. Both of these limitations are 
the result of the fact that selecting an origin of the coordinate system used to 
describe physical processes is always arbitrary. The nineteenth-century 
Scottish physicist James Clerk Maxwell (1831-1879) summarized the situation 
as follows.* 

We must, therefore, regard the energy of a material system as a quantity of 
which we may ascertain the increase or diminution as the system passes from 
one definite condition to another. The absolute value of the energy in the stan- 
dard condition is unknown to us, and it w r ould be of no value to us if we did 
know it, as all phenomena depend on the variations of energy and not on its ab- 
solute value. 

Next, we define the total energy of a particle to be the sum of the kinetic 
and potential energies: 

( 2 , 90 ) 

The total time derivative of E is 





*J. G. Maxwell, Matter and Motion (Cambridge, 1877), p. 91. 
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Substituting Equations 2.93 and 2.94 into 2.91, we find 



dE dU 

— = F-r + CVU) -f + — 
dt dt 



„ . dU 

— (F + VU) ’i H 

dt 



p = 



dU 

dt 



( 2 . 95 ) 



because the term F + V U vanishes in view of the definition of the potential en- 
ergy (Equation 2.88) if the total force is the conservative force F = —VC/. 

If C/is not an explicit function of the time (i.e., if dU/dt — 0; recall that we do 
not consider velocity-dependent potentials), the force field represented by F is 
conservative. Under these conditions, we have the third important conservation 
theorem: 



III. The total energy E of a particle in a consei'vative force field is a constant in time. 

It must be reiterated that we have not proved the conservation laws of linear 
momentum, angular momentum, and energy'. We have only derived various con- 
sequences of Newton’s laws; that is, if these laws are valid in a certain situation, 
then momentum and energy will be conserved. But we have become so enam- 
ored with these conservation theorems that we have elevated them to the-status 
of laws and we have come to insist that they be valid in any physical theory, even 
those that apply to situations in which Newtonian mechanics is not valid, as, for 
example, in the interaction of moving charges or in quantum-mechanical sys- 
tems. We do not actually have conservation laws in such situations, but rather 
conservation postulates that we force on the theory. For example, if we have two 
isolated moving electric charges, the electromagnetic forces between them are 
not conservative. We therefore endow the electromagnetic field with a certain 
amount of energy so that energy conservation will be valid. This procedure is sat- 
isfactory only if the consequences do not contradict any experimental fact, and 
this is indeed the case for moving charges. We therefore extend the usual con- 
cept of energy to include “electromagnetic energy” to satisfy our preconceived 
notion that energy must be conserved. This may seem an arbitrary’ and drastic 
step to take, but nothing, it is said, succeeds as does success, and these conserva- 
tion “laws” have been the most successful set of principles in physics. The refusal 
to relinquish energy and momentum conservation led Wolfgang Pauli 
(1900-1958) to postulate in 1930 the existence of the neutrino to account for 
the “missing” energy and momentum in radioactive (3 decay. This postulate al- 
lowed Enrico Fermi (1901-1954) to construct a successful theory of (3 decay in 
1934, but direct observation of the neutrino was not made until 1953 when 
Reines and Cowan performed their famous experiment.* By adhering to the 
conviction that energy and momentum must be conserved, a new elementary 



*C. L. Cowan, F. Reines, F. B. Harrison, H. W. Kruse, and A. D. McGuire, Science 124, 103 (1956). 
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particle was discovered, one that is of great importance in modern theories of 
nuclear and particle physics. This discovery is only one of the many advances in 
the understanding of the properties of matter that have resulted directly from 
the application of the conservation laws. 

We shall apply these conservation theorems to several physical situations in 
the remainder of this book, among them Rutherford scattering and planetary 
motion. A simple example here indicates the usefulness of the conservation 
theorems. 



EXAMPLE 2. 1 1 



A mouse of mass to J umps on the outside edge of a freely turning ceiling fan of 
rotational inertia / and radius R. By what ratio does the angular velocity 
change? 

Solution. Angular momentum must be conserved during the process. We are 
using the concept of rotational inertia learned in elementary physics to relate 
angular momentum L to angular velocity to: L = 1(0. The initial angular momen- 
tum L 0 = I(o 0 must be equal to the angular momentum L (fan plus mouse) 
after the mouse jumps on. The velocity of the outside edge is v, = wR. 

L = Ia> + mvR = — (J + mR l ) 

L = L 0 = Ioj ( j 

-(/+ to /? 2 ) = I- 

R y ' R 

v _ I 

v 0 I + mR 2 



w 0 / + to /? 2 



2.6 Energy 

The concept of energy was not nearly as popvdar in Newton’s time as it is today. 
Later we shall study two new formulations of dynamics, different from Newton’s, 
based on energy — the Lagrangian and Hamiltonian methods. 

Early in the nineteenth century, it became clear that heat was another form 
of energy and not a form of fluid (called “caloric”) that flowed between hot and 
cold bodies. Count Rumford* is generally given credit for realizing that the 



♦Benjamin Thompson (1753-1814) was born in Massachusetts and emigrated to Europe in 1776 as a 
loyalist refugee. Among the activities of his distinguished military and, later, scientific career, he su- 
pervised the boring of cannons as head of the Bavarian war department. 




2.6 ENERGY 



83 



great amount of heat generated during the boring of a cannon was caused by 
friction and not the caloric. If frictional energy is just heat energy, interchange- 
able with mechanical energy, then a total conservation of energy can occur. 

Throughout the nineteenth century, scientists performed experiments on 
the conservation of energy, resulting in the prominence given energy today. 
Hermann von Helmholtz (1821-1894) formulated the general law of conserva- 
tion of energy in 1847. He based his conclusion largely on the calorimetric ex- 
periments of James Prescott Joule (1818-1889) begun in 1840. 

Consider a point particle under the influence of a conservative force with 
potential U. The conservation of energy (actually, mechanical energy, to be pre- 
cise in this case) is reflected in Equation 2.90. 

E = T+ U = ~mv 2 + U(x) ( 2 . 96 ) 

where we consider only the one-dimensional case. We can rewrite Equation 2.96 
as 



v(t) 



dx 

dt 



~[E - V(x)] 
m 



( 2 . 97 ) 



and by integrating 



±dx 

U(x)] 



( 2 . 98 ) 



where x = x 0 at t = 1 0 . We have formally solved the one-dimensional case in 
Equation 2.98; that is, we have found x(t). All that remains is to insert the po- 
tential U(x) into Equation 2.98 and integrate, using computer techniques if 
necessary. We shall study later in some detail the potentials U — | kx 2 for har- 
monic oscillations and U ~ —k/xfor the gravitational force. 

We can learn a good deal about the motion of a particle simply by examin- 
ing a plot of an example of U(x) as shown in Figure 2-14. First, notice that, be- 
cause — T> 0, E> U(x) for any real physical motion. We see in Figure 2-14 
that the motion is bounded for energies £) and E 2 . For is), the motion is periodic 
between the turning points x a and x b . Similarly, for E 2 the motion is periodic, but 
there are two possible regions: x r s x :£ x d and x e ^ x ^ xy. The particle cannot 
‘jump” from one “pocket” to the other; once in a pocket, it must remain there 
forever if its energy remains at £ 2 . The motion for a particle with energy E () has 
only one value, x = x 0 . The pardcle is at rest with T = 0 [£ 0 = f/(x 0 ) ]. 

The motion for a particle with energy £ 3 is simple: The particle comes in 
from infinity, stops and turns at x = x g , and returns to infinity — much like a ten- 
nis ball bouncing against a practice wall. For the energy E 4 , the motion is un- 
bounded and the particle may be at any position. Its speed will change because it 
depends on the difference between E 4 and U(x). If it is moving to the right, it 
will speed up and slow down but continue to infinity. 
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U{x) 




FIGURE 2-14 Potential energy t/(x) curve with various energies E indicated. For 
certain energies, for example E x and £ 2 , the motion is bounded. 



The motion of a particle of energy E x is similar to that of a mass at the end 
of a spring. The potential in the region x a < x < x & can be approximated by 
U(x) = ~k{x - %o) 2 . A particle with energy barely above E 0 will oscillate about the 
point x = x 0 . We refer to such a point as an equilibrium point, because if the par- 
ticle is placed at x = x 0 it remains there. Equilibrium may be stable, unstable, or 
neutral. The equilibrium just discussed is stable because if the particle were 
placed on either side of x = x 0 it would eventually return there. We can use a 
hemispherical mixing bowl with a steel ball as an example. With the bowl right 
side up, the ball can roll around inside the bowl; but it will eventually settle to 
the bottom — in other words, there is a stable equilibrium. If we turn the bowl 
upside down and place the ball precisely outside at x = x 0 , the ball remains there 
in equilibrium. If we place the ball on either side of x = x 0 on the rounded sur- 
face, it rolls off; we call this unstable equilibrium. Neutral equilibrium would apply 
when the ball rolls on a flat, smooth, horizontal surface. 

In general, we can express the potential U(x) in a Taylor series about a 
certain equilibrium point. For mathematical simplicity, let us assume that the 
equilibrium point is at x = 0 rather than x = x 0 (if not, we can always redefine 
the coordinate system to make it so). Then we have 



lw “ lt + *UA + iiU r l + J, 



(2.99) 



The zero subscript indicates that the quantity is to be evaluated at x = 0. The po- 
tential energy U 0 at x = 0 is simply a constant that we can define to be zero without 
any loss of generality. If x = 0 is an equilibrium point, then 
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Equilibrium point 



and Equation 2.99 becomes 



U(x) 



x 2 fd*u\ x 3 / d $ U\ 

2! \dx 2 /o + 3! \dx 3 /o + 



( 2 . 100 ) 



( 2 . 101 ) 



Near the equilibrium point x = 0, the value of x is small, and each term in 
Equation 2.101 is considerably smaller than the previous one. Therefore, we 
keep only the first term in Equation 2.101: 





( 2 . 102 ) 



We can determine whether the equilibrium at x = 0 is stable or unstable by 
examining ( d 2 U/dx 2 ) 0 . If x = 0 is a stable equilibrium, U(x) must be greater 
(more positive) on either side of x = 0. Because x 2 is always positive, the con- 
ditions for the equilibrium are 





Stable equilibrium 
Unstable equilibrium 



( 2 . 103 ) 



If ( d 2 U/dx 2 ) 0 is zero, higher-order terms must be examined (see Problems 2-45 
and 2-46) . 



EXAMPLE 2.12 



Consider the system of pulleys, masses, and string shown in Figure 2-15. A light 
string of length b is attached at point A, passes over a pulley at point B located a 
distance 2d away, and finally attaches to mass m,. Another pulley with mass m 2 
attached passes over the string, pulling it down between A and B. Calculate the 
distance Xj when the system is in equilibrium, and determine whether the equi- 
librium is stable or unstable. The pulleys are massless. 

Solution. We can solve this example by either using forces (i.e., when 5q = 0 = 
Xj) or energy. We choose the energy method, because in equilibrium the ki- 
netic energy is zero and we need to deal only with the potential energy when 
Equation 2.100 applies. 

We let U = 0 along the line AB. 



U = — rri]gx } — + c) 



( 2 . 104 ) 
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B 2d A 




FIGURE 2- 1 5 Example 2.12. 



We assume that the pulley holding mass m 2 is small, so we can neglect the 
pulley radius. The distance cin Figure 2-15 is constant. 

x 2 = \/[(ft — x t ) 2 /4] — d~ 

U= —rriigx j — m.jg'V' [( b — x x ) 2 /^\ — d 2 — m^gc 

By setting dU/dXi — 0, we can determine the equilibrium position (x^o = x 0 : 

( dU\ , rn 2 g(b - x 0 ) 

= — rtiig H . = = 0 

\ dx Jo 4 V[(b - x„) 2 /4] - d 2 

4 m,y\^[(b — x 0 ) 2 /4] — d 2 = m 2 (6 — x 0 ) 

(6 — x 0 ) 2 (4 m\ — m^) — 16m\d 2 

4mi d 

*0 = b - 77== (2.10 

V 4wq — w 2 



Notice that a real solution exists only when 4mf > m|. 

Under what circumstances will the mass pull the mass mj up to the pul- 
ley JS(i.e., Xj = 0)? We can use Equation 2.103 to determine whether the equi- 
librium is stable or unstable: 

d 2 U ~m 2 g + m 2 g(b - x t ) 2 

dx \ 4{[(A - X074] - d 2 } m 16{[ (6 - X074] - d 2 f /2 



Now insert jq = x 0 . 



d 2 U \ _ g{ 4wf ~ ml) v 2 
dx f /,) 4 mid 



The condition for the equilibrium (real motion) previously was for 4 m? > mf, 
so the equilibrium, when it exists, will be stable, because (d 2 U/dx 2 ) 0 > 0. 
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EXAMPLE 2.18 



Consider the one-dimensional potential 



U(x) = 



- Wd 2 (x 2 + d 2 ) 
x 4 + 8rf 4 



( 2 . 106 ) 



Sketch the potential and discuss the motion at various values of x. Is the motion 
bounded or unbounded? Where are the equilibrium values? Are they stable or 
unstable? Find the turning points for E = — W/ 8. The value of Wis a positive 
constant. 

Solution, Rewrite the potential as 



Z(y) 



U{x) _ ~(y 2 + 1) 
W y 4 + 8 



where y = — 



( 2 . 107 ) 



First, find the equilibrium points, which will help guide us in sketching the 
potential. 

iZ -2; | 4;V + 1) 

dy y 4 + 8 (y 4 + 8) 2 

This is reduced to 



so 



y(y 4 + 2/ - 8) = 0 
y(y 2 + 4) (y 2 - 2) = 0 

yl = 2 , 0 

*bi ~ 0 

x 02 - V2 d > 

x 03 = — \fid t 



( 2 . 108 ) 



There are three equilibrium points. We sketch U(x)/ Wversus x/d in Figure 2-16. 
The equilibrium is stable at x 02 and x 08 but unstable at x 0l . The motion is 
bounded for all energies E < 0. We can determine turning points for any en- 
ergy E by setting E = U{x). 



E = 




- W(y 2 + 1 ) 
y 4 + 8 



( 2 . 109 ) 



y 4 + 8 = 8y 2 + 8 
y 4 - 8y 2 

y = ± 2V2, 0 ( 2 . 110 ) 

The turning points for E= — W/ 8 are x = — 2 V2rfand + 2 V2 d, as well as x — 0 — 
which is the unstable equilibrium point. 
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FIGURE 2-16 Example 2.13. Sketch of U{x)/W. 



2.7 Limitations of Newtonian Mechanics 

In this chapter, we have introduced such concepts as position, time, momentum, 
and energy. We have implied that these are all measurable quantities and that 
they can be specified with any desired accuracy, depending only on the degree 
of sophistication of our measuring instruments. Indeed, this implication appears 
to be verified by our experience with all macroscopic objects. At any given in- 
stant of time, for example, we can measure with great precision the position of, 
say, a planet in its orbit about the sun. A series of such measurements allows us to 
determine (also with great precision) the planet’s velocity at any given position. 

When we attempt to make precise measurements on microscopic objects, 
however, we find a fundamental limitation in the accuracy of the results. For 
example, we can conceivably measure the position of an electron by scattering a 
light photon from the electron. The wave character of the photon precludes an 
exact measurement, and we can determine the position of the electron only 
within some uncertainty Ax related to the extent (i.e., the wavelength) of the 
photon. By the very act of measurement, however, we have induced a change in 
the state of the electron, because the scattering of the photon imparts momen- 
tum to the electron. This momentum is uncertain by an amount A p. The product 
Ax A p is a measure of the precision with which we can simultaneously determine 
the electron’s position and momentum; Ax — > 0, A p — > 0 implies a measurement 
with all imaginable precision. It was shown by the German physicist Werner 
Heisenberg (1901-1976) in 1927 that this product must always be larger than a 
certain minimum value.* We cannot, then, simultaneously specify both the position 



*This result also applies to the measurement of energy at a particular time, in which case the prod- 
uct of the uncertainties is A E At (which has the same dimensions as Ax A p). 
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and momentum of the electron with infinite precision, for if Ax— »0, then we 
must have A p — » oo for Heisenberg’s uncertainty principle to be satisfied. 

The minimum value of Ax A p is of the order of 10 -34 J • s. This is extremely 
small by macroscopic standards, so for laboratory-scale objects there is no practi- 
cal difficulty in performing simultaneous measurements of position and momen- 
tum. Newton’s laws can therefore be applied as if position and momentum were 
precisely definable. But because of the uncertainty principle, Newtonian mechan- 
ics cannot be applied to microscopic systems. To overcome these fundamental 
difficulties in the Newtonian system, a new method of dealing with microscopic 
phenomena was developed, beginning in 1926. The work of Erwin Schrodinger 
(1887-1961), Heisenberg, Max Born (1872-1970), Paul Dirac (1902-1984), and 
others subsequently placed this new discipline on a firm foundation. Newtonian 
mechanics, then, is perfectly adequate for describing large-scale phenomena. But 
we need the new mechanics (quantum mechanics) to analyze processes in the 
atomic domain. As the size of the system increases, quantum mechanics goes over 
into the limiting form of Newtonian mechanics. 

In addition to the fundamental limitations of Newtonian mechanics as ap- 
plied to microscopic objects, there is another inherent difficulty in the 
Newtonian scheme — one that rests on the concept of time. In the Newtonian 
view, time is absolute, that is, it is supposed that it is always possible to determine 
unambiguously whether two events have occurred simultaneously or whether 
one has preceded the other. To decide on the time sequence of events, the two 
observers of the events must be in instantaneous communication,, either 
through some system of signals or by establishing two exactly synchronous clocks 
at the points of observation. But the setting of two clocks into exact synchronism 
requires the knowledge of the time of transit of a signal in one direction from one 
observer to the other. (We could accomplish this if we already had two synchro- 
nous clocks, but this is a circular argument.) When we actually measure signal 
velocities, however, we always obtain an average ve locity for propagation in oppo- 
site directions. And to devise an experiment to measure the velocity in only one 
direction inevitably leads to the introduction of some new assumption that we 
cannot verify before the experiment. 

We know that instantaneous communication by signaling is impossible: 
Interactions between material bodies propagate with finite velocity, and an inter- 
action of some sort must occur for a signal to be transmitted. The maximum ve- 
locity with which any signal can be propagated is that of light in free space: 
c = 3 x 10 8 m/s.* 

The difficulties in establishing a time scale between separate points lead us 
to believe that time is, after all, not absolute and that space and time are some- 
how intimately related. The solution to the dilemma was found during the pe- 
riod 1904-1905 by Hendrik Lorenz (1853-1928), Henri Poincare (1854-1912), 
and Albert Einstein (1879-1955) and is embodied in the special theory of rela- 
tivity (see Chapter 14). 



*The speed of light has now been defined to be 299,792,458.0 m/s to make comparisons of other 
measurements more standard. The meter is now defined as the distance traveled by light in a vac- 
uum during a time interval of 1 7299,792,458 of a second. 
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Newtonian mechanics is therefore subject to fundamental limitations when 
small distances or high velocities are encountered. Difficulties with Newtonian me- 
chanics may also occur when massive objects or enormous distances are involved. 
A practical limitation also occurs when the number of bodies constituting the sys- 
tem is large. In Chapter 8, we see that we cannot obtain a general solution in 
closed form for the motion of a system of more than two interacting bodies even 
for the relatively simple case of gravitational interaction. To calculate the motion in 
a three-body system, we must resort to a numerical approximation procedure. 
Although such a method is in principle capable of any desired accuracy, the labor 
involved is considerable. The motion in even more complex systems (for exam- 
ple, the system composed of all the major objects in the solar system) can like- 
wise be computed, but the procedure rapidly becomes too unwieldy to be of 
much use for any larger system. To calculate the motion of the individual mole- 
cules in, say, a cubic centimeter of gas containing ~ 10 19 molecules is clearly out 
of the question. A successful method of calculating the average properties of such 
systems was developed in the latter part of the nineteenth century by Boltzmann, 
Maxwell, Gibbs, Liouville, and others. These procedures allowed the dynamics 
of systems to be calculated from probability theory, and a statistical mechanics was 
evolved. Some comments regarding the formulation of statistical concepts in 
mechanics are found in Section 7.13. 



PROBLEMS 

2-1. Suppose that the force acting on a particle is factorable into one of the following 
forms: 

(a ) F(Xj, t) = f(x t )g(t) (b) F(Xi, t) = f(x t )g(t) (c) F(x„ x t ) = /(x, ■)£(*,■) 

For which cases are the equations of motion integrable? 

2-2. A particle of mass m is constrained to move on the surface of a sphere of radius R 
by an applied force F(0, <h ). Write the equation of motion. 

2-3. If a projectile is fired from the origin of the coordinate system with an initial veloc- 
ity t/ 0 and in a direction making an angle a with the horizontal, calculate the time 
required for the projectile to cross a line passing through the origin and making an 
angle /3 < a with the horizontal. 

2-4. A clown is juggling four balls simultaneously. Students use a video tape to deter- 
mine that it takes the clown 0.9 s to cycle each ball through his hands (including 
catching, transferring, and throwing) and to be ready to catch the next ball. What 
is the minimum vertical speed the clown must throw up each ball? 

2-5. A jet fighter pilot knows he is able to withstand an acceleration of 9g before black- 
ing out. The pilot points his plane vertically down while traveling at Mach 3 speed 
and intends to pull up in a circular maneuver before crashing into the ground, 
(a) Where does the maximum acceleration occur in the maneuver? (b) What is the 
minimum radius the pilot can take? 
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2-6. In the blizzard of ’88, a rancher was forced to drop hay bales from an airplane to 
feed her cattle. The plane flew horizontally at 160 kin/hr and dropped the bales 
from a height of 80 m above the flat range, (a) She wanted the bales of hay to land 
30 m behind the cattle so as to not hit them. Where should she push the bales out 
of the airplane? (b) To not hit the cattle, what is the largest time error she could 
make while pushing the bales out of the airplane? Ignore air resistance. 

2-7, Include air resistance for the bales of hay in the previous problem. A bale of hay 
has a mass of about 30 kg and an average area of about 0.2 m\ Let the resistance be 
proportional to the square of the speed and let c w = 0.8. Plot the trajectories with a 
computer if the hay bales land 30 m behind the catde for both including air resis- 
tance and not. If the bales of hay were released at the same time in the two cases, 
what is the distance between landing positions of the bales? 

2-8. A projectile is fired with a velocity u 0 such that it passes through two points both a 
distance h above the horizontal. Show that if the gun is adjusted for maximum 
range, the separation of the points is 




2-9. Consider a projectile fired vertically in a constant gravitational field. For the same 
initial velocities, compare the times required for the projectile to reach its maxi- 
mum height (a) for zero resisting force, (b) for a resisting force proportional to the 
instantaneous velocity of the projectile. 

2-10. Repeat Example 2.4 by performing a calculation using a computer to solve 
Equation 2.22. Use the following values: m — 1 kg, v 0 = 10 m/s, x a = 0, and k = 
0.1 s _1 . Make plots of ^versus t, x versus t, and v versus x. Compare with the results 
of Example 2.4 to see if your results are reasonable. 



2-11. Consider a particle of mass m whose motion starts from rest in a constant gravita- 
tional field. If a resisting force proportional to the square of the velocity (i.e., k mu' 2 ) 
is encountered, show that the distance s the particle falls in accelerating from n 0 to 
7.', is given by 



s(iy 0 ->-t/i) = — In 




2-12. A particle is projected vertically upward in a constant gravitational field with an 
initial speed v 0 . Show that if there is a retarding force proportional to the square 
of the instantaneous speed, the speed of the particle when it returns to the initial 
position is 



v o v t 




where v t is the terminal speed. 



2-13. A particle moves in a medium under the influence of a retarding force equal to 
mk{ v s + a 2 v), where k and a are constants. Show that for any value of the initial 
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speed the particle will never move a distance greater than rr/2ka and that the parti- 
cle comes to rest only for t-+ oo. 

2-14. A projectile is fired with initial speed v 0 at an elevation angle of a up a hill of slope 

P(a>i 8). 

(a) How far up the hill will the projectile land? 

(b) At what angle a will the range be a maximum? 

(c) What is the maximum range? 

2-15. A particle of mass m slides down an inclined plane under the influence of gravity. If 
the motion is resisted by a force / = kmv 2 , show that the time required to move a 
distance d after starting from rest is 

cosh - 1 (e ki ) 

%/ kg sin 0 

where 6 is the angle of inclination of the plane. 

2-16. A particle is projected with an initial velocity v 0 up a slope that makes an angle a 
with the horizontal. Assume frictionless motion and find the time required for the 
particle to return to its starting position. Find the time for v 0 = 2.4 m/s and a = 26°. 

2-17. A strong softball player smacks the ball at a height of 0.7 m above home plate. The 
ball leaves the player’s bat at an elevation angle of 35° and travels toward a fence 2 
m high and 60 m away in center field. What must the initial speed of the softball be 
to clear the center field fence? Ignore air resistance. 

2-18. Include air resistance proportional to the square of the ball’s speed in the previous 
problem. Let the drag coefficient be c w — 0.5, the softball radius be 5 cm and the 
mass be 200 g. (a) Find the initial speed of the softball needed now to clear the 
fence, (b) For this speed, find the initial elevation angle that allows the ball to most 
easily clear the fence. By how much does the ball now vertically clear the fence? 

2-19. If a projectile moves such that its distance from the point of projection is always in- 
creasing, find the maximum angle above the horizontal with which the particle 
could have been projected. (Assume no air resistance.) 

2-20. A gun fires a projectile of mass 10 kg of the type to which the curves of Figure 2-3 
apply. The muzzle velocity is 140 m/s. Through what angle must the barrel be ele- 
vated to hit a target on the same horizontal plane as the gun and 1000 m away? 
Compare the results with those for the case of no retardation. 

2-21. Show directly that the time rate of change of the angular momentum about the ori- 
gin for a projectile fired from the origin (constant g) is equal to the moment of 
force (or torque) about the origin. 

2-22. The motion of a charged particle in an electromagnetic field can be obtained from 
the Lorentz equation* for the force on a particle in such a field. If the electric field 
vector is E and the magnetic field vector is B, the force on a particle of mass m that 



*See, for example, Heald and Marion, Classical Electromagnetic Radiation (95, Section 1.7). 
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carries a charge q and has a velocity v is given by 

F = </E + qv X B 

where we assume that v <SC c(specd of light) . 

(a) If there is no electric field and if the particle enters the magnetic field in a di- 
rection perpendicular to the lines of magnetic flux, show that the trajectory is a 
circle with radius 



_ mv _ v 

qB (o c 

where <u f = qB/m is the cyclotron frequency. 

(b) Choose the z-axis to lie in the direction of B and let the plane containing E and 
B be the yz-plane. Thus 



B = Bk, E = E y j + E z k 

Show that the z component of the motion is given by 

qE. 

z(t) = z 0 + z 0 t + - -t 2 
2 m 



where 



z(0) = z 0 and z(0) = z 0 

(c) Continue the calculation and obtain expressions for x(t) and y(t ) . Show that the 
time averages of these velocity components are 

E y 

<*) = g- iy ) = 0 

(Show that the motion is periodic and then average over one complete period.) 

(d) Integrate the velocity equations found in (c) and show (with the initial condi- 
tions x(0) = -A/o) n x(0) = E y /B,y(0) = 0, y(0) = A) that 

_ A By A 

x{t) — cos u>J + — t, y{t) — — sin c o c t 

u) c B h) c 

These are the parametric equations of a trochoid. Sketch the projection of the 
trajectory on the ay-plane for the cases (i) A > \Ey/B\, (ii) A < \E y /B\, and 
(iii) A = | E y /B\ . (The last case yields a cycloid.) 



2-23. A particle of mass m = 1 kg is subjected to a one-dimensional force F(t) = kte~ al , 
where k — 1 N/s and a — 0.5 s~*. If the particle is initially at rest, calculate and plot 
with the aid of a computer the position, speed, and acceleration of the particle as a 
function of time. 



2-24. A skier weighing 90 kg starts from rest down a hill inclined at 17°. He skis 100 m 
down the hill and then coasts for 70 m along level snow undl he stops. Find the coef- 
ficient of kinedc fricdon between the skis and the snow. What velocity does the skier 
have at the bottom of the hill? 
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2-25. A block of mass m = 1.62 kg slides down a frictionless incline (Figure 2-A). The 
block is released a height h — 3.91 in above the bottom of the loop. 

(a) What is the force of the inclined track on the block at the bottom (point A)? 

(b) What is the force of the track on the block at point B ? 

(c) At what speed does the block leave the track? 

(d) How far away from point A does the block land on level ground? 

(e) Sketch the potential energy U(x) of the block. Indicate the total energy on the 
sketch. 



A 

h 



FIGURE 2-A Problem 2-25. 

2-26. A child slides a block of mass 2 kg along a slick kitchen floor. If the initial speed is 4 
m/s and the block hits a spring with spring constant 6 N/m, what is the maximum 
compression of the spring? What is the result if the block slides across 2 m of a 
rough floor that has fx k = 0.2? 

2-27. A rope having a total mass of 0.4 kg and total length 4 m has 0.6 m of the rope 
hanging vertically down off a work bench. How much work must be done to place 
all the rope on the bench? 

2-28. A superball of mass M and a marble of mass m are dropped from a height h with the 
marble just on top of the superball. A superball has a coefficient of restitution of 
nearly 1 (i.e., its collision is essentially elastic) . Ignore the sizes of the superball and 
marble. The superball collides with the floor, rebounds, and smacks the marble, 
which moves back up. How high does the marble go if all the motion is vertical? 
How high does the superball go? 

2-29. An automobile driver traveling down an 8% grade slams on his brakes and skids 30 
m before hitting a parked car. A lawyer hires an expert who measures the coeffi- 
cient of kinetic friction between the tires and road to be /x k = 0.45. Is the lawyer 
correct to accuse the driver of exceeding the 25-MPH speed limit? Explain. 

2-30. A student drops a water-filled balloon from the roof of the tallest building in town 
trying to hit her roommate on the ground (who is too quick). The first student 
ducks back but hears the water splash 4.021 s alter dropping the balloon. If the 
speed of sound is 331 m/s, find the height of the building, neglecting air resistance. 
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2-31. In Example 2,10, the initial velocity of the incoming charged particle had no com- 
ponent along the x-axis. Show that, even if it had an x component, the subsequent 
motion of the particle would be the same — that only the radius of the helix would 
be altered, 

2-32. Two blocks of unequal mass are connected by a string over a smooth pulley (Figure 
2-B) . If the coefficient of kinetic friction is fx what angle S of the incline allows the 
masses to move at a constant speed? 




FIGURE 2-B Problem 2-32. 



2-33. Perform a computer calculation for an object moving vertically in air under gravity 
and experiencing a retarding force proportional to the square of the object’s speed 
(see Equation 2.21). Use variables mfor mass and rfor the object’s radius; All the 
objects are dropped from rest from the top of a 100-m-tall building. Use a value of 
c K - — 0.5 and make computer plots of height y, speed v, and acceleration a versus t 
for the following conditions and answer the questions: 

(a) A baseball of m — 0.145 kg and r = 0.0366 m, 

(b) A ping-pong ball of m = 0.0024 kg and r = 0.019 m. 

(c) A raindrop of r — 0.003 m. 

(d) Do all the objects reach their terminal speeds? Discuss the values of the termi- 
nal velocities and explain their differences, 

(e) Why can a baseball be thrown farther than a ping-pong ball even though the 
baseball is so much more massive? 

(f) Discuss the terminal speeds of big and small raindrops. What are the terminal 
speeds of raindrops having radii 0.002 m and 0.004 m? 

2-34. A particle is released from rest (y = 0) and falls under the influence of gravity and 
air resistance. Find the relationship between v and the distance of falling y when 
the air resistance is equal to (a) av and (b) fiv 2 . 

2-35. Perform the numerical calculations of Example 2.7 for the values given in Figure 
2-8. Plot both Figures 2-8 and 2-9. Do not duplicate the solution in Appendix H; 
compose your own solution. 

2-36. A gun is located on a bluff of height h overlooking a river valley. If the muzzle ve- 
locity is v 0 , find the expression for the range as a function of the elevation angle of 
the gun. Solve numerically for the maximum range out into the valley for a given h 
and v fl . 
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2-37. A particle of mass m has speed v — a/x, where x is its displacement. Find the force 
F(x) responsible. 

2-38. The speed of a particle of mass m varies with the distance x as v{x) = ax - ". Assume 
i/(x = 0) = 0 at t = 0. (a) Find the force F(x) responsible, (b) Determine x(t) and 
(c) F(t). 

2-39. A boat with initial speed v 0 is launched on a lake. The boat is slowed by the water by 
a force F= —ae^ v . (a) Find an expression for the speed (b) Find the time and 
(c) distance for the boat to stop. 

2-40. A particle moves in a two-dimensional orbit defined by 

x(() = A(2at — sin at) 
y(t) = A(1 — cos at) 

(a) Find the tangential acceleration a t and normal acceleration a n as a function of 
time where the tangential and normal components are taken with respect to the 
velocity. 

(b) Determine at what times in the orbit a n has a maximum. 

2-41. A train moves along the tracks at a constant speed u. A woman on the train throws 
a ball of mass m straight ahead with a speed nwith respect to herself, (a) What is the 
kinetic energy gain of the ball as measured by a person on the train? (b) by a per- 
son standing by the railroad track? (c) How much work is done by the woman 
throwing he ball and (d) by the train? 

2-42. A solid cube of uniform density and sides of b is in equilibrium on top of a cylinder 
of radius R (Figure 2-C). The planes of four sides of the cube are parallel to the axis 
of the cylinder. The contact between cube and sphere is perfectly rough. Under 
what conditions is the equilibrium stable or not stable? 




FIGURE 2-C Problem 2-42. 

2-43. A particle is under the influence of a force F— —kx + kx^/a 2 , where k and a are 
constants and k is positive. Determine U(x) and discuss the motion. What happens 
when£= (l/4)ka 2 ? 

2-44. Solve Example 2.12 by using forces rather than energy. How can you determine 
whether the system equilibrium is stable or unstable? 
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2-45. Describe how to determine whether an equilibrium is stable or unstable when 
{d-U/dx 2 ) o = 0. 

2-46. Write the criteria for determining whether an equilibrium is stable or unstable 
when all derivatives up through order n, (d n U/dx n ) 0 = 0. 

2-47. Consider a particle moving in the region x > 0 under the influence of the potential 

, r I a x\ 

U{x) = f/ 0 - + - 
\x a) 

where U 0 — 1 J and a = 2 m. Plot the potential, find the equilibrium points, and 
determine whether they are maxima or minima. 

2-48. Two gravitationally bound stars with equal masses m , separated by a distance d, re- 
volve about their center of mass in circular orbits. Show that the period r is propor- 
tional to d 3/ 2 (Kepler’s Third Law) and find the proportionality constant. 

2-49. Two gravitationally bound stars with unequal masses mj and m 2 , separated by a dis- 
tance d, revolve about their center of mass in circular orbits. Show that the period r 
is proportional to d 3/2 (Kepler’s Third Law) and find the proportionality constant. 

2-50. According to special relativity, a particle of rest mass m 0 accelerated in one dimen- 
sion by a force F obeys the equation of motion dp/ dt — F. Here p — m 0 u/(l — 
v l /c 2 ) 1/2 is the relativistic momentum, which reduces to m 0 i> for xF/c 2 <?C 1. (a) For 
the case of constant /’and initial conditions x(0) = 0 = u(0), find x(f) and u(£). 
(b) Sketch your result for v(t). (c) Suppose that F/m 0 = 10 m/s 2 ( ~ g on Earth). 
How much time is required for the particle to reach half the speed of light and of 
99% the speed of light? 

2-51. Let us make the (unrealistic) assumption that a boat of mass m gliding with initial 
velocity v Q in water is slowed by a viscous retarding force of magnitude bv 2 , where b 
is a constant, (a) Find and sketch v{t). How r long does it take the boat to reach a 
speed of u o /1000? (b) Find x(t). How far does the boat travel in this time? Let m — 
200 kg, v 0 = 2 m/s, and b = 0.2 Nm _2 s 2 . 

2-52. A particle of mass m moving in one dimension has potential energy U{ x) — 
U 0 \_2(x/a ) 2 — (x/a) 4 ], where U n and a are positive constants, (a) Find the force 
F(x), which acts on the particle, (b) Sketch U{x). Find the positions of stable and 
unstable equilibrium, (c) What is the angular frequency u) of oscillations about the 
point of stable equilibrium? (d) What is the minimum speed the particle must have 
at the origin to escape to infinity? (e) At t = 0 the particle is at the origin and its ve- 
locity is positive and equal in magnitude to the escape speed of part (d). Find x(<) 
and sketch the result. 

2-53. Which of the following forces are conservative? If conservative, find the potential 
energy f/(r). (a) F x = ayz + bx + c, / v = ax z + bz, F z = axy + by. (b) F x — 

— ze~*, /’ v — lnz, F, = e~ x + y/z. (c) F = e r o/r(a, b, care constants). 

2-54. A potato of mass 0.5 kg moves under Earth’s gravity with an air resistive force of 

— kmv. (a) Find the terminal velocity if the potato is released from rest and k — 
0.01 s -1 . (b) Find the maximum height of the potato if it has the same value of k, 
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but it is initially shot directly upward with a student-made potato gun with an initial 
velocity of 1 20 m/s. 

2-55. A pumpkin of mass 5 kg shot out of a student-made cannon under air pressure at 
an elevation angle of 45° fell at a distance of 142 m from the cannon. The students 
used light beams and photocells to measure the initial velocity of 54 m/ s. If the air 
resistive force was F = —kmv, what was the value of k? 




CHAPTER 




Oscillations 



3.1 Introduction 

We begin by considering the oscillatory motion of a particle constrained to 
move in one dimension. We assume that a position of stable equilibrium exists 
for the particle, and we designate this point as the origin (see Section 2.6) . If the 
particle is displaced from the origin (in either direction), a certain force tends 
to restore the particle to its original position. An example is an atom in a long 
molecular chain. The restoring force is, in general, some complicated function 
of the displacement and perhaps of the particle’s velocity or even of some 
higher time derivative of the position coordinate. We consider here only cases 
in which the restoring force F is a function only of the displacement: F — F(x). 

We assume that the function F(x) that describes the restoring force possesses 
continuous derivatives of all orders so that the function can be expanded in a 
Taylor series: 



F(x) = F 0 + x 






(3.1) 



where F 0 is the value of F(x) at the origin (x = 0). and (d"F/dx”) 0 is the value of 
the nth derivative at the origin. Because the origin is defined to be the equilib- 
rium point, Tq must vanish, because otherwise the particle would move away from 
the equilibrium point and not return. If, then, we confine our attention to dis- 
placements of the particle that are sufficiently small, we can normally neglect all 
terms involving x 2 and higher powers of x. We have, therefore, the approximate 
relation 




(3.2) 
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where we have substituted k = — {di'/dx) () . Because the restoring force is always 
directed toward the equilibrium position (the origin), the derivative ( dF/dx) 0 is 
negative, and therefore k is a positive constant. Only the first power of the displace- 
ment occurs in F(x), so the restoring force in this approximation is a linear force. 

Physical systems described in terms of Equation 3.2 obey Hooke’s Law,* One 
of the classes of physical processes that can be treated by applying Hooke’s Law is 
that involving elastic deformations. As long as the displacements are small and the 
elastic limits are not exceeded, a linear restoring force can be used for problems 
of stretched springs, elastic springs, bending beams, and the like. But we must em- 
phasize that such calculations are only approximate, because essentially every real 
restoring force in nature is more complicated than the simple Hooke’s Law force. 
Linear forces are only useful approximations, and their validity is limited to cases in 
which the amplitudes of the oscillations are small (but see Problem 3-8) . 

Damped oscillations, usually resulting from friction, are almost always the 
type of oscillations that occur in nature. We learn in this chapter how to design an 
efficiently damped system. This damping of the oscillations may be counteracted 
if some mechanism supplies the system with energy from an external source at a 
rate equal to that absorbed by the damping medium. Motions of this type are 
called driven (or forced) oscillations. Normally sinusoidal, they have important 
applications in mechanical vibrations as well as in electrical systems. 

The extensive discussion of linear oscillatory systems is warranted by the 
great importance of oscillatory phenomena in many areas of physics and engi- 
neering. It is frequently permissible to use the linear approximation in the analy- 
sis of such systems. The usefulness of these analyses is due in large measure to 
the fact that we can usually use analytical methods. 

When we look more carefully at physical systems, we find that a large number 
of them are nonlinear in general. We will discuss nonlinear systems in Chapter 4. 



3.2 Simple Harmonic Oscillator 



The equation of motion for the simple harmonic oscillator may be obtained by 
substituting the Hooke’s Law force into the Newtonian equation F — ma. Thus 

— kx— mx ( 3 . 3 ) 



If we define 



wg = k/m 



( 3 . 4 ) 



Equation 3.3 becomes 

x + o>o* = 0 



( 3 . 5 ) 



* Robert Hooke (1635-1703). The equivalent of this force law was originally announced hy Hooke in 
1676 in the form of a Latin cryptogram: CEIIINOSSSTTUV. Hooke later provided a translation: ut 
tensio sic vis [the stretch is proportional to the force] . 
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According to the results of Appendix C, the solution of this equation can be 
expressed in either of the forms 

x(t) = A sin (w 0 l - 5) (3.6a) 

x(t) — A cos(d> 0 £ - <f>) (3.6b) 

where the phases* 8 and <f> differ by tt/2. (An alteration of the phase angle corre- 
sponds to a change of the instant that we designate t — 0, the origin of the time 
scale.) Equations 3.6a and b exhibit the well-known sinusoidal behavior of the 
displacement of the simple harmonic oscillator. 

We can obtain the relationship between the total energy of the oscillator 
and the amplitude of its motion as follows. Using Equation 3.6a for x{l), we find 
for the kinetic energy, 



T = 



- mx l = - mu> 2 A? cos 2 (w 0 t — 6) 
= ^ kA 2 cos 2 (a> Q t ~ 8) 



(3.7) 



The potential energy may be obtained by calculating the work required to 
displace the particle a distance x. The incremental amount of work dWnecessary 
to move the particle by an amount dx against the restoring force F is 

dW = —Fdx = kx dx . (3.8) 

Integrating from 0 to x and setting the work done on the particle equal to the 
potential energy, we have 



Then 




U — — kA 2 sin 2 (a»,)f — 8) 



(3.9) 



(3.10) 



Combining the expressions for T and U to find the total energy E, we have 
E= T+ U = ^kA 2 [cos 2 (oi 0 t - 8) + sin'Vot- 5)] 



E= T+ U = - kA 2 

2 



(3.11) 



so that the total energy is proportional to the squaw of the amplitude; this is a gen- 
eral result for linear systems. Notice also that E is independent of the time; that is, 



*The symbol d is often used to represent phase angle, and its value is either assigned or determined 
within the context of an application. Be careful when using equadons within this chapter because 5 
in one application may not be the same as the 5 in another. It might be prudent to assign subscripts, 
for example, S, and 5 2 , when using different equations. 
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energy is conserved. (Energy conservation is guar anteed, because we have been 
considering a system without frictional losses or other external forces.) 

The period r 0 of the motion is defined to be the time interval between succes- 
sive repetitions of the particle’s position and direction of motion. Such an inter- 
val occurs when the argument of the sine in Equation 3.6a increases by 2ir: 



or 



tt) 0 T 0 = 2 V 



(3.12) 




(3.13) 



From this expression, as well as from Equation 3.6, it should be clear that w 0 rep- 
resents the angular frequency of the motion, which is related to the frequency v 0 

by* 




(3.14) 



(3.15) 



Note that the period of the simple harmonic oscillator is independent of the 
amplitude (or total energy); a system exhibiting this property is said to be 

isochronous. 

For many problems, of which the simple pendulum is the best example, 
the equation of motion results in 8 + cofi sin d - 0, where 6 is the displacement 
angle from equilibrium, and o> 0 = Vg/€, where •£ is the length of the pendu- 
lum arm. We can make this differential equation describe simple harmonic 
motion by invoking the small oscillation assumption. If the oscillations about 
the equilibrium are small, we expand sin 0 and cos 0 in power series (see 
Appendix A) and keep only the lowest terms of importance. This often means 
sin 0 — Q and cos 0=1- 0 2 /2, where 0 is measured in radians. If we use the 
small oscillation approximation for the simple pendulum, the equation of mo- 
tion above becomes 0 + <wf0 = 0 , an equation that does represent simple har- 
monic motion. We shall often invoke this assumption throughout this text and 
in its problems. 



*1 lenceforlh we shall denote angular frequencies by o> (units: radians per unit time) and frequencies 
by v (units: vibrations per unit time or Hertz, Hz). Sometimes io will be referred io as a "frequency" 
for brevity, although “angular frequency” is to be understood. 





3.2 SIMPLE H ARMONIC OSCILLATOR 



103 



EXAMPLE 3.1 



Find the angular velocity and period of oscillation of a solid sphere of mass m 
and radius R about a point on its surface. See Figure 3-1. 



Solution. Let the rotational inertia of the sphere be l about the pivot point. In ele- 
mentary physics we learn that the value of the rotadonal inertia about an axis 
through the sphere s center is 2/5 mR 2 . If we use the parallel-axis theorem, the rota- 
tional inertia about the pivot point on the surface is 2/5 mR 2 + mR 2 — 7/5 mR 2 . 
The equilibrium position of the sphere occurs when the center of mass (center of 
sphere) is hanging directly below the pivot point. The gravitational force F = mg 
pulls the sphere back towards the equilibrium position as the sphere swings back 
and forth with angle 9. The torque on the sphere is N = la, where a = 9 is the an- 
gular acceleration. The torque is also N = R X F, with N— RF sin 6 = Rmg sin 6. 

For small oscillations, we have N = Rmgd. We must have 16 = —Rmg 6 for the 
equation of motion in this case, because as 6 increases, 9 is negative. We need to 
solve tire equation of motion for 9. 



6 + 




This equation is similar to Equation 3.5 and has solutions for the angular fre- 
quency and period from Equations 3,14 and 3.15, 




FIGURE 3-1 Example 3.1. The physical pendulum (sphere). 
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Note that the mass m does not enter. Only the distance R to the center of 
mass determines the oscillation frequency. 



3.3 Harmonic Oscillations in Two Dimensions 



We next consider the motion of a particle that is allowed two degrees of freedom. 
We take the restoring force to be proportional to the distance of the particle from 
a force center located at the origin and to be directed toward the origin: 



F = -kr 




(3.16) 


which can be resolved in polar coordinates into the components 




F x — — kr cos 8 = — kx 
F y — —kr sin 8 = ~ky t 


1 


(3-17) 


The equations of motion are 






X + WqX = Ol 
y + aily = ()J 




(3.18) 


where, as before, w 2 = k/rn. The solutions are 






x(<) = A cos (co n t — a)] 
y{t) — Bcos{w 0 t — (3) J 


■ 


(3.19) 



Thus, the motion is one of simple harmonic oscillation in each of the two direc- 
tions, both oscillations having the same frequency, but possibly differing in am- 
plitude and in phase. We can obtain the equation for the path of the particle by 
eliminating the time t between the two equations (Equation 3.19). First we write 

y{t) — B cos[w 0 ^ — a + (a — /3)] 

= B cos((o a l - a) cos(a — f3) — B sin(« n * — a)sin(a — (3) (3-20) 

Defining 8 = a — j3 and noting that cos(w 0 < — a) = x/A, we have 

y = ^ x cos 8 — B yjl — sin ® 



or 

Ay — Bx cos 5 — — ZiV 4 2 — x 2 sin 8 (3.21) 

On squaring, this becomes 

A 2 y' 2 — 2 ABxy cos 5 + B‘ z x 2 cos 2 8 = A' Z B 2 sin 2 8 - B z x 2 sin 2 8 

so that 

B 2 x 2 - 2 ABxy cos 8 + A 2 y 2 = A 2 B 2 sin 2 8 (3.22) 



3.3 HARMONIC OSCILLATIONS IN TWO DIMENSIONS 



If 5 is set equal to ± it/2, this equation reduces to the easily recognized equation 
for an ellipse: 



x 2 v 2 

- + -=! ; 8 = ±rr/2 



(3.23) 



If the amplitudes are equal, A= B, and if 8 = ±ir/ 2, we have the special case of 
circular motion: 



x 2 + y 2 = A' 1 , for A = B and 8 — ±v/2 (3.24) 

Another special case results if the phase 8 vanishes; then we have 
B 2 x 2 - 2ABxy + A 2 y 2 = 0 , 5-0 

Factoring, 

( Bx - Ay) 2 = 0 

which is the equation of a straight line: 

y = ~x, 5 = 0 (3.25) 

Similarly, the phase 5 = ±rr yields the straight line of opposite slope: 

B 

y= --X, 8= ±TT (3.26) 

A 

The curves of Figure 3-2 illustrate Equation 3.22 for the case A = B\8 — 90 ° 
or 270 ° yields a circle, and 5 = 180 ° or 360 °( 0 °) yields a straight line. All other 
values of 5 yield ellipses. 

In the general case of two-dimensional oscillations, the angular frequencies 
for the motions in the x- and y-directions need not be equal, so that Equation 
3.19 becomes 



(3.25) 



(3.26) 



x(t) = A cos((o x 
y(t) ~ B cos(a>j, 



J ~ «)1 

yt~P)l 



(3.27) 








S= 120° 



S = 150° 



S= 180° 



5 = 210° 








S = 240° 6 = 270° 8 = 300° 8 = 330° 8 = 360° 

FIGURE 3-2 Two-dimensional harmonic oscillation motion for various phase angles 
8 = a - p. 
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FIGURE 3-3 Closed two-dimensional oscillatory motion (called Lissajous curves) 
occurs under certain conditions for the xand y coordinates. 

The path of the motion is no longer an ellipse but a Lissajous curve.* Such a 
curve will be closed if the motion repeats itself at regular intervals of time. This 
will be possible only if the angular frequencies and w, are commensurable, that 
is, if <r) x /oj y is a rational fraction. Such a case is shown in Figure 3-3, in which 
(also a = /3). If the ratio of the angular frequencies is not a rational fraction, the 
curve will be open ; that is, the moving particle will never pass twice through the 
same point with the same velocity. In such a case, after a sufficiently long time 
has elapsed, the curve will pass arbitrarily close to any given point lying w’ithin 
the rectangle 2 A X 2B and will therefore "fill” the rectangle. r 

The two-dimensional oscillator is an example of a system in which an infini- 
tesimal change can result in a qualitatively different type of motion. The motion 
will be along a closed path if the two angular frequencies are commensurable. 
But if the angular frequency ratio deviates from a rational fraction by even an in- 
finitesimal amount, then the path will no longer be closed and it will “fill” the 
rectangle. For the path to be closed, the angular frequency ratio must be known to 
be a rational fraction with infinite precision. 

If the angular frequencies for the motions in the x- and y-directions are dif- 
ferent, the shape of the resulting Lissajous curve- strongly depends on the phase 
difference 8 = a — (3. Figure 3-4 shows the results for the case <w = 2<o x tor 
phase differences of 0, tt/3, and 7t/2. 



3.4 Phase Diagrams 

The state of motion of a one-dimensional oscillator, such as that discussed in 
Section 3.2, will be completely specified as a function of time if two quantities 



♦The French physicist Jules Lissajous (1822-1880) demonstrated this in 1857 and is generally given 
credit, although Nathaniel Bowditch seems to have reported in 1815 two mutually orthogonal oscil- 
lations displaying the same motion (Cr81 ). 

+A proof is given, for example, by Haag (Ha62, p. 86). 
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ti\ = 261,. ; 8 = n/3 




co y = 2 o ) x ; S = it / 2 



FIGURE 3-4 Lissajous curves depend strongly on the phase differences of the angle 5. 



are given at one instant of time, that is, the initial conditions x(t 0 ) and x(/ 0 ). 
(Two quantities are needed because the differential equation for the motion is of 
second order.) We may consider the quantities x(l) and x(() to be the coordinates of 
a point in a two-dimensional space, called phase space. (In two dimensions, the 
phase space is a phase plane. But for a general oscillator with n degrees of freedom, 
the phase space is a 2«-dimensional space.) As the lime varies, the point P(x, x) 
describing the state of the oscillating particle will move along a certain phase path in 
the phase plane. For different initial conditions of the oscillator, the motion will be 
described by different phase paths. Any given path represents the complete time his- 
tory of the oscillator for a certain set of initial conditions. The totality of all possible 
phase paths constitutes the phase portrait or the phase diagram of the oscillator.* 
According to the results of the preceding section, we have, for the simple har- 
monic oscillator, 

x(l) — A sin(tti 0 / — 8) (3.28a) 

x(l) = Aa>(i cos (to n l — 8) (3.28b) 



If we eliminate t from these equations, we find for the equation of the path 





(3.29) 



This equation represents a family of ellipses,* several of which are shown in Figure 3-5. 
We know' that the total energy £ of the oscillator is (Equation 3.11), and be- 
cause <>)% = k/m. Equation 3.29 can be written as 



2 E/k + 2 E/m 



= 1 



(3.30) 



Each phase path, then, corresponds to a definite total energy' of the oscillator. This 
result is expected because the system is conservative (i.e., £ = const.). 

No two phase paths of the oscillator can cross. If they could cross, this would 
imply that for a given set of initial conditions x(t 0 ), x(t Q ) (i.e., the coordinates of the 



*These considerations are not restricted to oscillating particles or oscillating systems. The concept of 
phase space is applied extensively in various fields of physics, particularly in statistical mechanics. 
fThc ordinate of the phase plane is sometimes chosen to be x/<o {) instead of x; the phase paths are 
then circles. 
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FIGURE 3-5 Phase diagram for a simple harmonic oscillator for a variety of total 
energies E. 



crossing point), the motion could proceed along different phase paths. But this is 
impossible because the solution of the differential equation is unique. 

If the coordinate axes of the phase plane are chosen as in Figure 3-5, the 
motion of the representative point P(x, x) will always be in a clockwise direction, 
because for x > 0 the velocity x is always decreasing and for x < 0 the velocity is 
always increasing. 

To obtain Equations 3.28 for x(t) and x(t ) , we must integrate Equation 3.5, a 
second-order differential equation: 

d^x 

~^ + 0>lx= 0 (3.31) 



We can obtain the equation for the phase path, however, by a simpler procedure, 
because Equation 3.31 can be replaced by the pair of equations 



dx 




dx 

Jt = 



(3.32) 



If we divide the second of these equations by the first, we obtain 



dx 

dx 



~o)l 



x 

X 



(3.33) 



This is a first-order differential equation for x = x(x ) , the solution to which is just 
Equation 3.29. For the simple harmonic oscillator, there is no difficulty in obtaining 
the general solution for the motion by solving the second-order equation. But in more 
complicated situations, it is sometimes considerably easier to directly find the equation 
of the phase path x = x(x) without proceeding through the calculation of x(t) . 



3.5 Damped Oscillations 

The motion represented by the simple harmonic oscillator is termed a free oscilla- 
tion; once set into oscillation, the motion would never cease. This oversimplifies 
the actual physical case, in which dissipative or frictional forces would eventually 
damp the motion to the point that the oscillations would no longer occur. We can 
analyze the motion in such a case by incorporating into the differential equation a 
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term representing the damping force. It does not seem reasonable that the damping 
force should, in general, depend on the displacement, but it could be a function 
of the velocity or perhaps of some higher time derivative of the displacement. It is 
frequently assumed that the damping force is a linear function of the velocity,* 
¥ d = av. We consider here only one-dimensional damped oscillations so that we 
can represent the damping term by — bx. The parameter b must be positive in order 
that the force indeed be resisting. (A force —bx with b < 0 would act to increase the 
speed instead of decreasing it as any resisting force must.) Thus, if a particle of 
mass m moves under the combined influence of a linear restoring force — kx and a 
resisting force — bx, the differential equation describing the motion is 



?nx + bx + kx — 0 



(3.34) 



w'hich we can write as 



x + 2/3x + (o^x = 0 



(3.35) 



Here (3 = b/2m is the damping parameter and to f) = \fk/m is the characteristic 
angular frequency in the absence of damping, The roots of the auxiliary equation 
are (cf. Equation C.8, Appendix C) 

r { = -jS + Vf3- - «iij | 

/ f (3-36) 

r 2 = - j8 - \//3 2 - wjj J 

The general solution of Equation 3.35 is therefore 



x(t) = e exp(V/3 2 — to^t) + A 2 exp( — V/3 2 — ;)] 



(3.37) 



There are three general cases of interest: 

Underdamping: wjj > /3 2 

Critical damping: a>l = /3 2 

Overdamping: a>l < 

The motion of the three cases is shown schematically in Figure 3-6 for specific initial 
conditions. We shall see that only the case of underdamping results in oscillatory 
motion. These three cases are discussed separately. 

Underdamped Motion 

For the case of underdamped motion, it is convenient to define 

io\ = u>l - p 1 (3.38) 



See Section 2.4 for a discussion of the dependence of resisting forces on velocity. 
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FIGURE 3-6 Damped oscillator motion for three cases of damping. 



where <wf > 0; then the exponents in the brackets of Equation 3.37 are imaginary, 
and the solution becomes 

x{t) = <T0'[ ( 3 . 39 ) 
Equation 3.39 can be rewritten as* 

x(t) = Ae~0 ( cos(co l t — S ) ( 3 . 40 ) 

We call the quantity oj t the angular frequency of the damped oscillator. Strictly • 
speaking, we cannot define a frequency when damping is present, because the 
motion is not periodic — that is, the oscillator never passes twice through a given 
point with the same velocity. However, because co } — 27t/( 27’ ] ), where 1\ is the 
time between adjacent zero x-axis crossings, the angular frequency w, has meaning 
for a given time period. Note that 27j would be the “period" in this case, not 7j. 
For simplicity, we refer to aq as the “angular frequency” of the damped oscillator, 
and we note that this quantity is less than the frequency of the oscillator in the ab- 
sence of damping (i.e., w, < tt> 0 ). If the damping is small, then 

= Vo>0 - j3 2 = <O 0 

so the term angular frequency may be used. But the meaning is not precise unless 

0 = 0 . 

The maximum amplitude of the motion of the damped oscillator decreases 
with time because of the factor exp( — /3 1), where f3 > 0, and the envelope of the 
displacement versus time curve is given by 

x en = ±Ae~f >< ( 3 . 41 ) 

This envelope and the displacement curve are shown in Figure 3-7 for the case 5 = 0. 
The sinusoidal curve for undamped motion (j3 = 0) is also shown in this figure. A 
close comparison of the two curves indicates that the frequency for the damped 
case is less (i.e., that the period is longer) than that for the undamped case. 



*Sec Exercise l)-6, Appendix D. 
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FIGURE 3-7 The underdamped motion (solid line) is an oscillatory motion (short 
dashes) that decreases within the exponential envelope (long dashes). 



The ratio of the amplitudes of the oscillation at two successive maxima is 



Ae-W 
Ae~* T+ T '> 



e^i 



(3.42) 



where the first of any pair of maxima occurs at : t = T and where t, = 2tt/<» x . The 
quantity exp(/3T,) is called die decrement of die modon; the logarithm of exp (/Ur , ) — 
that is, /3r, — is known as the logarithmic decrement of the motion. 

Unlike die simple harmonic oscillator discussed previously, the energy of the 
damped oscillator is not constant in time; rather, energy is continually given up to 
the damping medium and dissipated as heat (or, perhaps, as radiadon in the form 
of fluid waves). The rate of energy loss is propordonal to the square of the velocity 
(see Problem 3-11), so the decrease of energy docs not take place uniformly. The 
loss rate will be a maximum w r hen the parucle attains its maximum velocity near 
(but not exactly at) the equilibrium position, and it will instantaneously vanish 
when the particle is at maximum amplitude and has zero velocity. Figure 3-8 shows 
the total energy and the rate of energy loss for the damped oscillator. 



EXAMPLE 3.2 



Construct a general phase diagram analytically for the damped oscillator. Then, 
using a computer, make a plot for x and x versus t and a phase diagram for the 
following values: A = 1 cm, o> 0 = 1 rad/s, ft — 0.2 s" 1 , and 8 — v/2 rad. 

Solution. First, we write the expressions for the displacement and the velocity: 

x(t) = Ae ^cos(te)[f— 8) 

x(t) = ~Ae~ e> [fi cos (<U]f — 5) + <W] sin(co 1 f — 5)] 

These equations can be coverted into a more easily recognized form by introducing 
a change of variables according to the following linear transformations: 

u = u)\X , w = (3x + x 



u — ( 0 [Ae cosiw^ — 8) 
w = sin(w 1 ? — 8) 



Then 
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FIGURE 3-8 The total energy and rate of energy loss for the damped oscillator. 




FIGURE 3-9 Example 3.2. 



If we represent u and aiin polar coordinates (Figure 3-9), then 



P 



\/ ii 1 + u>' 2 , 



4> = d>-\t 



Thus 



p = <i i x Ae~^ /u ^ 



which is the equation of a logarithmic spiral. Because the transformation from *, 
k to u, w is linear, the phase path has basically the same shape in the u-w plane 
(Figure 3-1 Oa) and x-x plane (Figure 3-10b) . They both show a spiral phase path 
of the underdamped oscillator. The continually decreasing magnitude of the radius 
vector for a representative point in the phase plane always indicates damped 
motion of the oscillator. 
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x,x 




Time (s) 

(c) 

FIGURE 3-10 Results for Example 3.2. The phase path (a) of the W y u coordinates 
and (b) of the x, x coordinates, and (c) a numerical calculation of 
position and speed versus time. The spiral path is characteristic of the 
underdamped oscillator. 



The actual calculation using numbers can be done by various means with a 
computer. We chose to use one of the commercially available numerical pro- 
grams that has good graphics output. We chose the values A = 1, /3 = 0.2, h = 1, 
m = 1, and 8 = tt/2 in the appropriate units to produce Figure 3-10. For the 
particular value of 8 chosen, the amplitude has x= 0 at t = 0, but x has a large 
positive value, which causes x to rise to a maximum value of about 0.7 m at 2 s 
(Figure 3-10c). The weak damping parameter (i allows the system to oscillate 
about zero several times (Figure 3-1 0c) before the system finally spirals down to 
zero. The system crosses the x = 0 line eleven times before x decreases finally to 
less than 10 -3 of its maximum amplitude. The phase diagram of Figure 3-10b 
displays the actual path. 






